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ABSTRACT 

We present a detailed imaging and spectroscopic study of the conjugate hard X-ray (HXR) footpoints 
(FPs) observed with RHESSIm the 2003 October 29 XIO flare. The double FPs first move toward and 
then away from each other, mainly parallel and perpendicular to the magnetic neutral line, respectively. 
The transition of these two phases of FP unshearing motions coincides with the direction reversal of 
the motion of the loop-top (LT) source, and with the minima of the estimated loop length and LT 
height. The FPs show temporal correlations between HXR flux, spectral index, and magnetic fleld 
strength. The HXR flux exponentially correlates with the magnetic fleld strength, which also anti- 
correlates with the spectral index before the second HXR peak's maximum, suggesting that particle 
acceleration sensitively depends on the magnetic fleld strength and/or reconnection rate. Asymmetries 
are observed between the FPs: on average, the eastern FP is 2.2 times brighter in HXR flux and 1.8 
times weaker in magnetic fleld strength, and moves 2.8 times faster away from the neutral line than 
the western FP; the estimated coronal column density to the eastern FP from the LT source is 1.7 
times smaller. The two FPs have marginally different spectral indexes. The eastern-to- western FP 
HXR flux ratio and magnetic fleld strength ratio are anti-correlated only before the second HXR 
peak's maximum. Neither magnetic mirroring nor column density alone can explain the totality of 
these observations, but their combination, together with other transport effects, might provide a full 
explanation. We have also developed novel techniques to remove particle contamination from HXR 
counts and to estimate effects of pulse pileup in imaging spectroscopy, which can be applied to other 
RHESSI flares in similar circumstances. 

Subject headings: acceleration of particles — Sun: flares — Sun: magnetic flelds — Sun: X-rays 



1. INTRODUCTION 

In the c l assica l picture o f solar flares (CHSKP, 
Carmichaell _ ll964 [Hiravamal Il974l : ISturrockl Il966l : 
Kopp fc Pneuman Il976l ). magnetic reconnection occurs 
high in the corona, resulting in plasma heating and 
particle acceleration. Some of the high-energy par- 
ticles may be trapped in the corona by plasma tur- 
bulence (e.g., iPetrosian fc Liul l2004f) an d /or b y ma^- 
netic mirroring (e.g.. iMelrose fc White! 1 19791 : iLeachI 
[1984), where energetic electrons pro duce the hard X-ray 
(HXR) l oop-top (LT) sources (e.g., "Masuda et al."l994'; 
[Petrosia n et all [20 02: Sui fc Holman 2003; Jiang et al. 
[2006; B attadia fc Benzl 120061 : iLiul [2006). Some par- 
ticles escape the acceleration region into interplane- 
tary space and can contribute to solar ener g:etic par- 
ticle (SEP) events detect ed at 1 AU (e.g., iLiu et all 
l2004al : iKrucker et~ani2007[ ). A signiflcant portion of the 
accelerated particles escape down magnetic fleld lines 
and produce bremsstrahlung HXRs (by electrons) or 
gamma-ray lines ( by protons and other ions), primarily 
in the thick-target ([Br ownl [19711 : [Petrosianlll973t) transi- 
tion region and chromosphere. This results in the com- 
monly observed HXRs at t he con jugate footpoints (FPs; 
iHovng et al. 1981: Saka ol Il994l : [Petrosian et al. 2002; 
ISaint-Hilaire et al.ll2008l ) of the flare loop. Consequent 
energy redistribution in the lower atmosphere along the 
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arcade of such loops produces ribbons seen in Ha and 
occasionally in white-light. As time proceeds, reconnec- 
tion occurs at higher altitudes in the inverted-Y s haped 
confi guration, a nd consequently the HXR FPs ([Sakaol 
"1994': > Qiu etal] [20021 : [Liu et al.ll2004b[ ) and Ha ribbons 
(^Svestka[[l976[ ) are usually seen to move away from each 
other, while the soft and hard X-ray LT sources are ob- 
served to move upward simultaneo usly at a comparable 
spee d (Gahagher et al. 2002; Sui et~al] [20041 : [Liu et al.[ 
[2004bl ). At the same tin ie, an upward g:rQwin gs: loop sys- 
tem can be seen in Ha irm k Tajiakal [19735 . extreme 
ultravi olet (EUV, [G allaghe r et al.[[2002[). and soft X-ray 
fSXR: IPallavicini et al.Hl975[ : [Tsuneta et al.[[T992[ ). De- 
tails of reconnection and particle acceleration, however, 
remain largely unknown. X-ray observations of the LT 
and FP sources, particularly of their spatial, temporal, 
and spectral properties, combined with magnetic fleld 
measurements of the flare region, can provide critical in- 
formation about how and where electrons are accelerated 
subsequent to magnetic reconnection. 

Apparent motions of X-ray LT and FP sources can be 
understood as sequential excitat ions /formations of flare 
loops (e.g.. [Schrijver et al] [2006 1) when the primary re- 
connection site changes its location. Source motions 
deviating from the above classical picture have been 
observed, especially in the past decade: (1) The alti- 
tude of the HXR LT source was discovered by the Ra- 
maty High Energy Solar Spectroscopic Imager {RHESSI) 
to decrease during the rise of the impulsive phase 
prior to the expected increase (ISui Ho lman 2Q0l; 
[Sui et al.[ [20041 : [Liu et al.[ [2004aL [2008[ ). Shrinkage of 
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loops was also observed in SXR (Svestka et al.l [TqSTI: 
Forbes fc Actonl Il996l: [Vrsnak et al l l2QQ6l: iReeves et al ' 
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and microwave (|Li fc Gan 



2QQ5|) during flar es, and in SXR on the quiet Sun 
(jWang et al. 1997). ( 2) HXR FPs exhibit even more pat- 
terns of motion. Sa kao et al.l (p. 998) found that in 7 out 
of 14 flares observed by the Yohkoh Hard X-Ray Tele- 
scope (HXT) the FPs move away from each other, while 
in the rest of the flares the FP separation decreases (see 
also Fletcher fc Hudson"2002) or remains roughly con- 
stant. Bogachev et al. (2005) extended this study to 31 
Yohkoh flares and found that only 13% (type I) are con- 
sistent with the classical flare model showing double FPs 
moving away from and nearly perpendicular to the mag- 
netic neutral line (NL), while 26% (type II) of the flares 
exhibit FPs moving mainly along the NL in antiparal- 
lel directions, 35% (type III) have F Ps moving along the 
NL in the same dir ection (see also iKrucker et al.l 120031 : 
iGrigis fc Benj[2Q05f ). and the remaining 26% show com- 
plex motion patterns. Their type II events are of par- 
ticular interest as they suggest that flare loops excited 
or formed more recently during reconnection are less 
sheared^ (e.^., see a revie w on Hinitori results bv lTanakal 
119871 : iSomov et all l2002[ ). Reduction of shear during 
flares has been observed for decades in the form of loops 
seen in Ha (Zirin fc Tanaka 1973; Rust fc Bar 1973) and 
EUV (jAsai et al.ll2003 ) with an increasing angle from the 
NL, or of app arent unshear ing motions of FPs seen in 
Ra (lA^aLeLa h 2003), EUV (iSu et al.l[2QQ7h. and HXR s 
(jSakaol 11994: Masuda eTaT "2001"; "Schmahl et alj [2QQ6h . 
(3) Recently Ji et al. (2006, 2007) found that the de- 
crease of the (Ha and HXR) FP separation and shear 
occurred when the estimated height of the LT source de- 
creases during the rise of the impulsive peak. This has 
spurred renewed interest by establishing the connection 
between the LT descending motion and the FP approach- 
ing motion with decreasing shear. LT descents reported 
in the past were usually observed in flares occurring near 
the solar limb where the LT height can be readily mea- 
sured, but the FP motions in the east-west direction are 
obscured by projection effects. Flares close to disk cen- 
ter, like the one reported here, give an alternative per- 
spective. 

Correlations between a pair of conjugate HXR FPs 
are expected, since they are believed to be produced by 
high-energy electrons released from the same accelera- 
tion region. The relative timing of conjugate FPs was 
found to be simult aneous within an uncertainty of 0.1- 
0.3 s (|Sakaol 119941 ) based on Yohkoh Hard X-Ray Tele- 
scope (HXT) observations. For double FPs in tens of 
flares observed by the Ramaty High Energy Solar Spec- 
troscopic Imager {RHESSI), temporal correlations in the 
HXR fluxes in two wide energy bands (25-50 and 50- 
10 keV) with a tim e resolution of 8 s were identifled 
bv I Jin fc Dind (|2007l ). Spectral correlations at individ- 



^ Hagyard et al.l (| 19841 ) defined the degree of magnetic shear as 
"the difference at the photosphere between the azimuths of a po- 
tential field and the observed vector field, where the potential field 
is the one satisfying the boundary conditions provided by the ob- 
served line-of-sight field" . In this paper we use "shear" to mean the 
angle between the line connecting the two conjugate FPs (i.e., be- 
lieved to be on the same magnetic loop) and the normal to the NL; 
our "shear angle" defined in this sense represents apparent shear, 
since we do not have accurate vector magnetic field measurements. 



ualHXR peaks were investigated by ISaint-Hilaire et aP 
(j2008[ ). who found power-law indexes that differed by 
<0.6. This spec tral in dex difference is similar to that 
found by Saka^ (|1994| ), b ut smaller than t he va lues as 
high as 1 or 2 reported bv iPetrosian et al] (|2002l ). both 
based on analysis of Yohkoh HXT images. 

Asymmetric FPs, i.e., conjugate FPs with differ- 
ent properties (HXR fluxes, magn etic fleld str engths, 
etc.) are commonly observed (e.g.. [Sakaol 11994 ). This 
has been ascribed to asymmetric magnetic mirroring 
where a brighte r HXR FP is associated with a weaker 
magnetic fleld ([Li et al. 1 19971: Aschwanden et al.l 1 19991 : 



iQiu et all [2MI : iLi fc Dind [2QQ1 . This picture is con- 
sistent with observations at radio wavelengths where 
brighter microwave emission appe ars at the FP with 
the stronger ma gnetic fleld (e.g., iKundu et al.l 1 19951 : 
IWang et al.lll995l). Exceptions to the mirroring scenario 
were reported bv lGoff et al.l (|2004l ), who found one third 
of 32 Yohkoh flares with an opposite trend, that is, the as- 
sociation of the brighter HXR FP wi t h the stronger mag- 
netic fleld. Falewicz fc Siarkowskil (|2007[ ) re-examined 
three exceptions in the sample of Goff et al. and at- 
tributed this opposite asymmetry to different column 
densities i n the two leg s of the flare loo p, as also sug- 
gested by Emslie et aO (|2003[ ) and (|2006D . Tem- 
poral variatio ns of the flux asymmetry we re found in a 
Yohkoh flare (|Siarkowski fc Falewicj [20041 ) . and energy- 
and t ime-dependent var iations were seen in a RHESSI 
flare (|Alexander fc Met calf 20Q 2|). Th e latter were in- 
terpreted by "McClements fc Alexander (2005) as a con- 
sequence of an asymmetric, energy- and time-dependent 
injection of accelerated electrons. 

Previous studies of conjugate HXR FPs, in general, 
suffered from limited time, spatial, and/or energy resolu- 
tion and/or coverage of HXR emission, mainly restricted 
by the instrumental capabilities, or from lack of magnetic 
fleld measurements. We report here on a comprehen- 
sive study of the conjugate FPs in the 2003 October 29 
XI flare observed by RHESSI that overcomes many of 
the previous shortcomings. This flare provides a unique 
opportunity to track the spatial and spectral evolution 
of the double HXR FPs and their associated magnetic 
flelds in great detail, and to study all three interrelated 
aspects: unshearing motions^ correlations^ and asymme- 
tries. This flare occurred near disk center, where FP mo- 
tions and line-of-sight magnetic fleld measurements have 
minimum projection effects. Its long (~20 minutes) im- 
pulsive phase and high RHESSI count rates up to several 
hundred keV allow for a detailed study of variations both 
in time and energy. The flare was also well observed by 
the Transition Region and Coronal Explorer {TRACE) ^ 
the Solar and Heliospheric Observatory {SOHO)^ and 
other spacecraft and many ground-based observatories. 
The rich database of multiwavelength observations and 
a wide range o f literature cov e ring different aspec ts of 
this event fe.g.. IXu et alll2004l : [Krucker et al.ll2005f ) are 
particularly beneflcial for this in-depth study. 

We present the observations and data analysis in § O 
These include general RHESSI light curves, images, and 
imaging spectroscopy. We investigate in § [3] the two 
phases (fast and slow) of unshearing motions of the FPs 
and the associated LT motion. In § [H we explore var- 
ious correlations of the FPs, particularly of their HXR 
fluxes, spectral shapes, spatial variations, and magnetic 
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fields. Possible contributions to the HXR flux and spec- 
tral asymmetries are discussed in § [3 followed by our 
conclusions in § [6l A discussion of pulse pileup effects, 
technical details on coaligning images made by different 
instruments, a mathematical treatment of the asymmet- 
ric column density effect, and an estimate of the coronal 
column densities in the legs of the flare loop are given 
in Appendixes [Al [Bl [Cl and [Dl respectively. Note that 
due to the particular circumstances of this flare, we em- 
ployed new procedures to remove particle contamination 
from HXR light curves and to estimate effects of pulse 
pileup in imaging spectroscopy. These techniques, as a 
bonus of this paper, can be applied to other RHESSI 
flares in similar difficult circumstances that could not be 
analyzed before. 

2. OBSERVATIONS AND DATA ANALYSIS 

We present in this section general multiwavelength and 
RHESSI X-ray observations to indicate the context for 
our detailed discussions to follow on the conjugate FPs. 
The event under study occurred in AR 10486 (W5°S18°) 
starting at 20:37 UT on 2003 O ctober 29, during 
the so-called Halloween storms (e.g.. ' Gopalswamy et al.l 
12005). It was a Geostationary Operational Environ- 
ment Satellite {GOES) XIO class, white-light, two- 
ribbon flare, whi ch pro duced strong gamma-ray line 
emission (Hurfor d et al.l [2006) and helioseismic signals 
(|Donea fc Lindsevll2005f ). It was associated with various 
other solar activity, including a fast (~2000 kms~^) halo 
coronal mass ejections (CMEs), and heliospheric conse- 
quences. This was the first white-light flare observed at 
the opacity minimum at 1.6 /im, which corresponds to 
the deepest layer of the photosphere that can be seen 
(|Xu et al.l[20M . I2006f ). There were strong photospheric 
shearing flows present near the magnetic NLs in this ac - 
tive region prior to the flare onset ( Yan g~et al.l 12004 ). 
which may be related to the unusually lar ge amount o f 
free magnetic energy (~6 x 10^^ ergs; M etcalf et al.l[2005h 
stored in this A R. By analyzing Huairou and Mees vector 
magnetograms. lLiu et al.l (|2007 ) proposed that this flare 
resulted from reconnection between magnetic flux tubes 
having opposite current helicities. This may be con- 
nected to the soft X-ray sigmoid structure and unshear- 
ing motions of HXR FPs foun d by[jLet all (2008) during 
the early phase of the flare. Liu fc Havashii (^006), us- 
ing potential field extrapolations from the SOHO Michel- 
son Doppler Imager (MDI) observations, investigated the 
large-scale coronal magneti c field of AR 104 86 and its 
high productivity of CMEs. iLiu et all (|2006ar ) found re- 
mote brightenings more than 2x10^ km away from the 
main flare site. Solar energetic particles (SEPs) were de- 
tected after this flare by GOES and the Advanced Com- 
position Explorer (ACE). 

Our goal in this paper is to understand the tempo- 
ral and spectral variations of the asymmetric HXR FPs 
and their associated magnetic fields. We thus focus on 
HXR observations obtained by RHESSI and line-of-sight 
(LOS) photospheric magnetograms obtained by SOHO 
MDI. Vector magnetograms measured with chromo- 
spheric emission lines are more desirable for this study, 
as relevant magnetic mirroring may take place above the 
chromosphere where thick-target HXRs are produced. 
The only available vector magnetograms recorded dur- 
ing/near the flare time are Mees data, which, however. 



has not yet been fully calibrated (K. D. Leka 2008, pri- 
vate communication) and thus is not used here. On 
the other hand, since this flare is close to disk center 
(W5°S18°), LOS MDI magnetograms provide a reason- 
able approximation of the photospheric field that is as- 
sumed to be proportional to the chromospheric field. Yet 
bear in mind that the field is likely tilted and can ob- 
scure the MDI measurement due to projection effects. 
It would have been interesting to examine microwave 
images which may sh ow opposite FP asymmetry as in 
HXRs (|Kundu et al.lfT995 ). However, spectrograms of 
this fla re obtained at the Ow ens Valley Solar Array 
(OVSA iGarv fc HurfoTdl flOOOh do not allow for image 
reconstruction due to poor data quality (J. Lee fc C. Liu 
2007, private communication), while Nobeyama was not 
observing (before local time 6 AM). 

2.1. Light Curves and X-ray Images 

RHESSI had very good coverage of this event. 
However, HXR counts, particularly at high energies 
(>40 keV), were contaminated by counts produced by ra- 
diation belt particles bombarding the spacecraft from all 
directions. Particle-caused counts mainly influence the 
rear segments of RHESSI detectors (unshielded), while 
the front segments are less affected during the flare since 
they have a smaller volume (^0.14 that of the rear seg- 
ments) and they stop most of the down ward- incident 
flare photons <150 keV (|Smith et al.ll200 2^ that may ex- 
ceed the number of particle counts. In addition, parti- 
cle count rates are not modulated by the grids and ap- 
pear as a DC offset in the modulation pattern, which 
is rem oved during image reconstruction ([Hurford et al.l 
I2OO2I ). Therefore, particles do not affect our analysis in 
this paper since it mainly relies on images made using 
front rate only. 

We obtained the spatially integrated X-ray fluxes free 
from particle contamination in the following way. (1) 
At low energies (<30 keV), no correction is needed since 
ample fiare counts dominate over particle counts. (2) At 
high energies (30-150 keV), we used the rear segment 
counts to estimate the contribution of particles to the 
front segments, since rear counts in this energy range are 
almost entirely produced by particles. To do this, we 
first selected a non-flare interval 20:42:40-20:47:40 UT 
on 2003 October 30, one day after the flare when the 
spacecraft was at approximately the same geomagnetic 
location, to get a background estimate. We obtained the 
background count rate spectra averaged over detectors 
1, 3-6, 8, and 9 for front and rear segments separately, 
and calculated the front-to-rear count-rate ratio. (The 
ratio is close to the volume ratio ~0.14 of the two seg- 
ments, because particle bombardments are expected to 
be isotropic. Meanwhile, the weak energy dependence 
of the ratio, ranging from 0.08 at 30 keV to 0.14 at 150 
keV, may be related to the geometry of the segments.) 
We then repeated this for the flare to accumulate front 
and rear segment count rate spectra for every 4 s interval 
from 20:37 to 21:07 UT. For each time and energy bin, 
the rear count rate was multiplied by the front-to-rear 
ratio at this energy obtained above and then subtracted 
from the front count rate. A sample of the count rates 
before and after this correction is shown in Figure [H As 
expected, during the bulk of the flare duration between 
20:40 and 21:00 UT, the estimated fractional particle 
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contamination is minimal (<15% of total counts, except 
up to 50% between 20:52 and 20:55 UT) at a lower energy 
(55-56 keV), and becomes more appreciable (up to 75% 
at -20:53:30 UT) at a higher energy (120-125 keV). Note 
that this technique cannot be used for energies <30 keV 
due to the threshold (lower-level discrimi nator) of the 
rear segments set at -20 keV (|Smith et al.ll2002f ). 




20:40 20:45 20:50 20:55 21:00 21:05 
Time (since 2003/10/29 20:37) 

Fig. 1. — Count rates at two selected energies averaged over 
front segments 1, 3-6, 8, and 9 before (dotted) and after (solid) the 
correction for particle contamination. The 120-125 keV curves are 
scaled by a factor of 1/200. 

The resulting X-ray count rates are also integrated 
in wide energy bands and are shown in Figure [2] to- 
gether with GOES soft X-ray and OVSA microwave 
fluxes. We flnd that the RHESSI 80-120 keV and OVSA 
16.4 GHz fluxes follow each other closely in time before 
20:55 UT and both exhibit two major peaks (Peaks 1 
and 2) divided at 20:48 UT. Note that the frequency 
at the maximum of the OVSA microwave spectrum is 
<11.2 GHz (except for three 4 s intervals between 20:41 
and 20:44 UT when it reaches up to 14 GHz), and thus 
16.4 GHz is on the optically-thin side of the spectrum 
which can be fltted with a power law. 



TABLE 1 

Phases during the course of the flare. 





Phase I 


Phase n 




20:40:40-20:44 20:44-20:59:40 


FP unshearing motion 


fast 


slow 


FP motion w.r.t. NL 


parallel 


perpendicular 


LT altitude (estimated) 


decrease 


increase 




Peak 1 


Peak 2 




20:39-20:48 


20:48-21:02 


mag. mirroring asymmetry strong 


weak 


column densities in loops 


small 


large 


FP B-fields correlation 


exists 


disappears 



To obtain the flare morphology and its general evo- 
lution, we focused on a time range from 20:40:40 to 
20:59:40 UT beyond which the double conjugate FPs of 
interest (identifled below) were not clearly imaged, due 
to complex morphology and/or low count rates. We flrst 
divided this time range into 57 consecutive 20 s inter- 
vals, except for one interval that was shortened to 12 s 
to avoid the decimation state change at 20:46:36 UT. We 
then reconstructed images in two broad energy bands, 
12-25 and 60-100 keV, using the CLEAN algorithm and 
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HXR Peak 1 ■ HXR Peak 2 




20:40 20:45 20:50 20:55 21:00 21:05 
Time (since 2003/10/29 20:37:00) 



Fig. 2. X-ray fluxes from RHESSI (6-12 and 80-120 keV, left 
scale) and GOES-12 (1-8 A, right scale), and microwave flux (16.4 
GHz, in sfu, left scale) from OVSA. The RHESSI Huyies are average 
count rates of front segments 1, 3-6, 8, and 9, with contamination 
from radiation belt particles removed. The artiflcial steps are due 
to attenuator state changes. The OVSA 16.4 GHz and RHESSI 
80-120 keV curves are scaled by factors of 0.06 and y, respectively. 
The vertical dotted lines mark the two transitions, ti and t2, where 
ti divides Phases I and H based on FP motions and t2 divides two 
major HXR peaks (see Table Count rates in the 20-50 keV 
range are affected pulse pileup and thus are not shown here. 

unifo rm weighting among detectors 3-8 (|Hurford et alJ 
I2002D . The effective FWHM angular resolution is 9.8''. 

A sample of the resulting images is shown in Figure [3l 
Early in the flare (before 20:43:20 UT, Fig.[3H), several 
bright points at 60-100 keV are dispersed across the im- 
age, suggesting FPs of multiple loops. Part of the 12- 
25 keV emission appears elongated and curved between 
the adjacent FPs, corresponding to the LT source(s). 
Toward the south-west, part of the 12-25 keV emission 
seems to overlap with the FP emission, possibly due to 
a projection effect. As time proceeds, the FP structure 
seen at 60-100 keV becomes simpler, and only two dis- 
tinct FPs are present (after 20:43:20 UT). They generally 
move away from each other. At the same time, the 12- 
25 keV emission gradually changes from one to two LT 
sources, one in the north and the other in the south. 

We identified the conjugate FPs and the correspond- 
ing LT source of interest as follows for detailed anal- 
ysis: (1) At later times (after 20:43:20 UT), only two 
FPs are seen in each image at 60-100 keV and so they 
are considered conjugate. We call the FP on the east- 
ern (left) side E-FP and the one on the western (right) 
side W-FP (see, e.g.. Fig. [Spi). (2) At earlier times when 
more than two FPs are present, we set forth the fol- 
lowing selection criteria: (a) The source morphology of 
the two conjugate FPs must be consistent with the pic- 
ture that they are magnetically connected through the 
LT source between them seen in the corresponding 12- 
25 keV image (see, e.g.. Fig. (b) During the time 

evolution the two FPs must show continuity and con- 
sistency in position and HXR flux, which other short- 
lived FPs lack. Under these criteria, the selected E-FP 
is the brightest FP to the east of the magnetic NL {thick 
dashed in Figs. ^ and[3a), and W-FP is the one to the 
west located nearest to the NL. (3) Once the conjugate 
FPs are found, their corresponding LT source was iden- 
tified as the 12-25 keV emission that lies closest to the 
straight line joining the FPs. For example, at later times 
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□ [2005/10/29 20:40:40-20:41:00] 
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Fig. 3. — Evolution of HXR sources as seen in negative RHESSI CLEAN images made with detectors 3-8 (FWHM resolution 9.8^0 
in two energy bands, 12-25 keV as gray scale and 60-100 keV as contours (at 23%, 35%, 60%, and 90% of the maximum of each image). 
The dashed line in panel h shows the simplified magnetic NL also plotted in Fig. [Jli. Correction for the solar differential rotation has been 
applied to this NL, and source centroids and images presented elsewhere in this paper when applicable unless otherwise noted. 



(see, e.g., Fig.[3Fi), the northern LT is selected, while the 
southern LT is ignored since it does not seem to have any 
corresponding FP emission, presumably because of its 
faintness t hat exceeds RHESSFs dynam ic range (>10:1 
for images. [Hurford et al.l[2QQ2l : lLiull2QQ6l . p. 214 therein). 

2.2. Imaging Spectroscopy of Footpoint and Loop-top 

Sources 

Next, we examine the spectroscopic characteristics of 
the LT and FP sources and their temporal evolution. 
For each of the 57 consecutive 20 s intervals defined 
above, we reconstructed CLEAN images in 16 energy 
bins that are progressively wider from 6 to 150 keV. A 
sample of these images is shown in Figure H] for 20:51:20- 
20:51:40 UT, where four images showing similar morphol- 
ogy as in neighboring energy bins are omitted. The emis- 
sion is dominated by the two LT sources at low energies 
and the double FP sources at high energies. 

The next step was to obtain photon fluxes of the 
sources for each time interval. For each FP source, we 
used a hand-drawn polygon that envelops all the 10% 
(of the maximum brightness of the image) contours at 
energies where this FP source was clearly imaged. For 
the corresponding LT source, we drew a polygon that 
encloses the 20% contours, which was selected to mini- 
mize spatial contamination from the FPs. We then read 
the resulting multiple-energy image cube into the stan- 
dard RHESSI spectral analysis software (OSPEX) pack- 
age. This package integrates photon fluxes inside each 
polygon, and uses the full detector response matrix to 
estimate the true incident photon spectrum. The RMS 
of the residual map of the CLEAN image was used to cal- 
culate the uncertainty for the photon flux in each energy 
bin, with proper consideration of the source area and grid 



spatial resolu tion. This imaging spectroscopy technique 
is detailed bv lLiu et al.l (2008). Note that we did not use 
contours at a fixed level (as opposed to polygons fixed in 
space) to obtain the fluxes because of the complex source 
morphology that makes such contours vary with energy. 

One important i ssue for this XIO flare is pulse pileup 
(|Smith et al.ll2002l ) that at high count rates distorts the 
count-rate spectrum. We have discussed in Appendix [A] 
various effects of pileup on our analysis and the rem- 
edy that we have applied to minimize them. Although 
it is currently not possible to obtain accurate spectra 
throughout the full energy range for all sources, pileup 
mainly affects the LT spectra in the range of 20-50 keV 
(e.g., see Figs.Spi and[4^c). In other words, pileup effects 
on spectral shapes are negligible for the LT sources be- 
low 20 keV and for the FP sources above 50 keV. This 
conclusion enabled us to confine the extent of pileup ef- 
fects both in energy and in space. We thus fitted the 
LT spectrum below 20 keV wi th an assumed isot hermal 
model from CHIANTI ver. 5.2 (jYoung et al.ll2003[ ). using 
the default coronal iron abundance of 4 times the pho- 
tospheric value, to determine its temperature (T) and 
emission measure (EM); we fitted the FP spectrum above 
50 keV with an assumed single power-law model to find 
its spectral index (7) and normalization flux (/) at the 
reference energy of 50 keV. 

Spectra of the LT and FP sources are shown in Fig- 
ure [5]a for the interval of 20:44:40-20:45:00 UT (during 
the main impulsive peak). Above 50 keV, both FP spec- 
tra have a power-law shape, with the E-FP flux being 
twice that of W-FP but only slightly harder. Conse- 
quently, the W-to-E ratio of the two FP spectra generally 
decreases with energy (Fig. ^) or stays constant within 
uncertainties. Below 20 keV the LT spectrum shows the 
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Fig. 4. — negative CLEAN images in different energy bins at 20:51:20-20:51:40 UT made with detectors 3-9. The contour levels are 
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cm~-^ s"-*^ arcsec"^), of each individual image. The number below /max indicates the total counts accumulated by the detectors used. The 
same three dark, hand-drawn polygons in each panel were used to obtain the fluxes of the LT and two FP sources. 
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FlG. 5. — (a) Spatially resolved spectra of the LT and two 
FP sources at 20:44:40-20:45:00 UT. The horizontal bars on the 
LT spectrum represent the energy bin widths in the range of 6- 
150 keV. The dashed and dotted lines are power-law flts from 50- 
150 keV for E-FP and W-FP, respectively, and the dot-dashed line 
is a single temperature thermal flt from 6-19 keV for the LT. The 
legend shows the photon fluxes (/) at 50 keV and the spectral 
indexes (7) for the FPs, and the emission measure (EM) and tem- 
perature (T) for the LT. (6) Ratio of the W-to-E FP fluxes, (c) 
Fitting residuals normalized by the la uncertainties for the LT 
(dot- dashed), E-FP (dashed) and W-FP (dotted) sources. 

exponential shape of isothermal bremsstrahlung emis- 
sion, with the iron line feature at 6.7 keV visible. Note 
that below 50 keV the FP spectra may be compromised 
by pileup effects^ and spatial contamination from the LT 

^ We note that the same power-law trends of the two FP spectra 
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Fig. 6. — Evolution of spectroscopic parameters of the LT source. 
(a) Temperature T and (b) emission measure EM with 20 s inte- 
gration time. The solid lines are 7-point box-car smooths of the 
original fltting results indicated by the s ymb ols, (c) Radius r of the 
equivalent sphere of the LT source (see ^2.2|) and (d) corresponding 
electron number density tilt- 

source, and likewise above 20 keV the apparent LT flux 
is contaminated by FP emission at the same energy and 
by pileup from lower energies (Fig.[5]a). We show in Fig- 
ure [6] the spectral evolution of the LT source and defer 
that of the FP sources to § KT\ 
In order to infer the density of the LT source, we as- 

extend below 50 keV down to ~30 keV, suggesting that pulse pileup 
may have minimal effects on the spectral shapes of the FPs, and 
that our selection of 50 keV as the lower limit for reliable FP spectra 
is likely to be unnecessarily conservative. 
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sumed that it has a spherical shape with the projected 
area a equal to the area inside the 50% brightness contour 
at 12-25 keV. We then obtained the radius r = (a/7r)^/^ 
and volume V = Anr^/S of the equivalent sphere and the 
corresponding density tilt = [EM/(V/)]^/^. In doing 
so we assumed a filling factor / of unity, which means 
the density obtained here would be a lower limit, and 
used the EM values smoothed with a 7-point box-car to 
minimize fluctuations possibly caused by the inevitable 
anti-correlation between T and EM during spectral fit- 
ting. The values of r and tilt as functions of time are 
shown in Figures [6]c and[6H, respectively. As evident, the 
size of the sphere stays roughly constant between 11''- 
16'' and thus the density follows the same trend as the 
EM. 

3. TWO-PHASE FOOTPOINT UNSHEARING AND 
LOOP-TOP MOTIONS 

We now examine in detail the spatial evolution of the 
double FP sources and the corresponding LT source, by 
tracking the migration of their emission centroids. For 
each 12-25 image obtained in § 12.11 we used a contour 
at 50% of the maximum brightness of the LT source to 
locate its centroid, while for each 60-100 keV image we 
used a 90% contour of each conjugate FP. The reason for 
a higher contour level for the FPs (than the LT) is that 
the E-FP source spreads along the flare ribbon (see, e.g.. 
Fig. [3]e) and we need this brightest "kernel" to obtain 
the corresponding magnetic field strength at the FP (see 
^4.3p . The resulting centroids are shown in Figure [3 

The background preflare MDI magnetogram in Fig- 
ure [Th was corrected from SO HO LI view to Earth view 
and shifted by Ax = 4.5" ± 2.0" and Ay = -2.8" ± 2.0" 
in the solar east-west (x) and south-north (y) direc- 
tions, respectively, to match the RHES SI dispect believed 
to be have sub-arcsecond accuracy (Fi vian et al. 2002). 
The required shifts were determined b y cross-correlating 
MDI magnetic anom aly features (e.g. JOiu fc Garvll2003l : 
iSchriiver et"alll2006[ ) with HXR FPs, as described in Ap- 
pendix [Bl The MDI map and all RHES SI centroids were 
corrected for the solar rotation and shifted to their cor- 
responding positions at a common time (20:50:42 UT) in 
the middle of the flare. As is evident, E-FP is located in 
the negative {dark) polarity to the left of the simplified 
magnetic NL {red dashed)^ while W-FP is in the positive 
{white) polarity to the right. 

In an attempt to correct for projection effects and to 
obtain the true 3D loop geometry, we assumed that the 
centroids of the LT and two FP sources at a given time 
are connected by a semi-circular loop. We then used the 
solar X and y coordinates of these three points in the 
sky plane to determine the size and the orientation of 
the semi-circle in 3D space, knowing that the FPs are 
located on the solar surface and the LT in the corona. 
A sample of the loops at selected times is shown in Fig- 
ure [3). We find that the inclination angle between the 
model loop and the vertical direction ranges from 14° to 
63°, and that the loop length (/total ; see Fig. [Sla) gener- 
ally first decreases and then increases with a minimum 
at 20:43:50 UT. 

The LT centroids {plus signs) as shown in Figure [3) 
are situated at all times close to the NL {red dashed) as 
expected, and form two clusters, one in the south and the 
other in the north. As time proceeds, the LT centroid ap- 



pears to move from the apex of one loop to another along 
the arcade seen in TRACE 195 A (not shown). It first 
gradually moves southward until 20:43:30 UT (marked 
by the middle circle in Fig.[7l?)), when it starts to rapidly 
shift to the northern cluster and then continue moving 
northward at progressively lower speeds. This can be 
more clearly seen from its relative displacement projected 
onto the north-south NL as a function of time (Fig. [8]?), 
AyLr)' The height of the LT centroid estimated from 
the model loops is also shown in Figure [8]?) and exhibits 
a general decrease followed by an increase. 

As to the FPs, in general, E-FP first moves south- 
ward and then turns to the east, while W-FP first moves 
northward and then turns to the west, as indicated by 
the thick arrows in Figure [3). The evolution of the 
position of E-FP relative to W-FP is shown in Fig- 
ure [Tic. There is clearly a turning point which occurs at 
ti =20:44 UT and divides the evolution of the FP posi- 
tions into two phases: (1) Phase I (20:40:40-20:44:00 UT) 
when the two FPs generally move toward each other in 
a direction essentially parallel to the NL, (2) Phase II 
(20:44:00-20:59:40 UT) when the FPs move away from 
each other main l y perp endicular to the NL. According to 
iBogachev et al.l (|2005l ). this fiare falls into their type II 
events during Phase I and then type III during Phase II. 
Another signature of this two-phase division is the mor- 
phological transition at 20:43:20 UT, before which there 
are multiple FP sources, but only two FPs present after- 
wards (see Fig. [3]). Below we describe in detail the HXR 
source evolution in the two phases. 

We further decomposed the distance between the FPs 
into two components: perpendicular and parallel to the 
NL as shown in Figure Elc, where the two phases are 
divided by the vertical dotted line at ti. As can be seen, 
the parallel distance {asterisks) first rapidly decreases at 
a velocity of 108 ±18 kms~^ given by the linear fit during 
Phase I; it then stays almost constant during Phase II 
with a slow increase (7 ± 1 kms~^). In contrast, the 
perpendicular distance {squares) has a slow variation in 
Phase I (15 ± 13 kms~^) and increases continuously at 
a velocity of 51 ± 1 kms~^ in Phase II. These velocities 
are comparable to those of the TRACE EUV FPs found 
by ISchrijver et"al.i (^06) in the 2003 October 28 X17 
flare that occurred in the same active region as the flare 
under study. From the speed (~60 km s~^) and size 
(<1400 km) of EUV FPs, they inferred the crossing time 
of the FP diameter or the excitation time scale of the 
HXR-producing electron beam in a single flare loop to 
be <23 s. 

Next we obtained the shear angle {6; Fig. [71c) from the 
normal to the NL (parallel to the ?/-axis) to the W-to-E 
relative positional vector, which is shown as a function 
of time in Figure [8](i. This angle exhibits a fast decrease 
(from 56° to 22°) during Phase I and a slow decrease 
(down to 12°) during Phase II. An independent study by 
iJi et al.l (|2008l ). with different identifications of the FPs 
in this flare and thus larger scatter, also found a simi- 
lar decrease of the shear angle in two phases, which they 
referred to as sigmoid and arcade phases based on the X- 
ray morphology. The apparent unshearing motions of the 
HXR FPs indicate that the later reconnected magnetic 
field lines are less sheared. It can be seen that TRACE 
195 A loops corresponding to the HXR FPs at early 
times (not shown) are indeed highly sheared. Similar 



8 



Liu et al. 






20 40 60 


80 

X (orcsecs 


100 120 
) 


140 160 


-360 






1 neutral line 




-570 






\ 1 




u -380 


i-(b) Source centroids: 








-390 


60-100 keV: E-FP A 
^ W-FPv 
12-25 keV: LT + 




1\ °- 


cl 


-400 


; 20:42 20:46 20:50 20:54 


20:58 






40 60 


80 

X (orcsecs] 


100 


120 




-60 -40 
AX (orcsecs) 



W-FP 



Fig. 7. — (a) preflare SOHO MDI magnetogram (centered at 
20:30:35 UT with an integration time of 30 s) overlaid with mag- 
netic NLs in white and the centroids of RHESSI FPs (E- and W- 
FP, at 60-100 keV) in color. The white (dark) gray scale represents 
positive (negative) line-of-sight magnetic fields pointing away (to- 
ward) the observer. The vertical red dashed line here and in (6) 
is the simplified NL that follows the observed magnetic NL be- 
tween the two FPs. The temporal evolution at 20 s intervals of the 
RHESSI centroids is indicated by the color bar in (6). (6) Evolu- 
tion of the centroid of the RHESSI LT source at 12-25 keV {plus 
signs). The tracks of the FP centroids shown in (a) are repeated 
here as dotted lines, along with positions at four selected times 
(labeled a, d, /, and i in Fig. [3]) marked by triangles. We show 
semi-circular model loops at these times projected onto the sky 
plane, each of which connects three centroids (of two FPs & one 
LT) of the same color. The thick dark arrows indicate the general 
direction of motion for the LT and two FPs. (c) Relative centroid 
positions of E-FP with respect to W-FP which is selected as the 
origin. The start and end of the time evolution are marked by the 
black and red arrows. is the shear angle between the normal 
(due west) to the NL and the line joining the two FP centroids. 
The open circle here and those in (6) mark the transitional time 

=20:44 UT from fast to slow unshearing motions of the FPs (see 

unsh earing motions w ere observed i n various wavelengths 
fe^ lZirin fc Tanakal ll973: Ma suda et alJl2QQll : [Siret all 
I2QQ7I ) . Note that an opposite process took place prior to 
the flare, that is, strong photospheric shearing flows ob- 
served near the NL (Yang et al. 2004, as mentioned ear- 
lier). This process increased the shear of fleld line s and 
built up magnetic stress and free energy ( Metcalf et al.l 
l2005f ) in the system during the preflare phase. 




20:44 20: 
Time 



48 20:52 20:56 21:00 

(UT, 2003/10/29) 



Fig. 8. — (a) History of distance (/i and I2) from the LT cen- 
troid to the centroids of the two FPs along the semi-circular model 
loop as shown in Fig. [Tt, together with the length of the loop 
(^totai = h +^2)- {b) Relative displacement (A^/lt, dot-dashed) of 
the LT centroid parallel to the magnetic NL shown in Fig. and 
estimated height of the LT centroid above the solar surface (solid). 
(c) Orthogonal components of the separation between W-FP and 
E-FP perpendicular (squares) and paral\e\(asterisks) to the NL. 
The straight lines are linear fits to the distances in the two phases, 
labeled with the corresponding velocities (in kms"-*^). (d) Shear 
angle 6 defined in Fig. [Tic. (e) Perpendicular distances to the NL 
from E-FP (diamonds) and W-FP (crosses). Linear fits for the 
whole fiare duration are shown as the straight lines. The error bar 
shown on the last data point of each line here and in (c) is the 
RMS deviation of the data from the corresponding fit. 



Finally we investigate the relationship between the FP 
and LT motions. The division at ti =20:44 UT between 
the two phases of the FP unshearing motions coincides 
(within 30 s) with the minimum of the estimated loop 
length (Fig. [8]a) and the direction reversal of the ap- 
parent LT motion (Fig. ^) noted above. Also the es- 
timated LT height undergoes a general decrease during 
Phase I when the HXR flux is on the rise. A similar 
decrease of the LT altitude during the rising portion 
of the impulsive phase, followed by a subsequent alti- 
tude increase, has been observed in many RHESSI B^diies 
near the l imb (e . g. ISui fc Holmanll2003l: ISmelall l2004l : 
iLiu et all l2004al . l2008l : iHolman et all l2005f ) . For those 
flares a complete physical picture is obscured because the 
observed FP motions are strongly subject to projection 
effects, but this drawback vanishes for the disk flare un- 
der study here. Assuming that our semi-circular model 
loops yield reasonable estimates of the LT heights and 
loop lengths, the above source motions, when taken to- 
gether, suggest the following scenario: (1) During Phase 
I, as the the reconnection site and thus the LT source 
migrate southward along the NL or the arcade, shorter 
and less sheared loops are energized, which leads to the 
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apparent decrease of the LT altitude and the fast un- 
shearing motion of the FPs. (2) During Phase II, as the 
the reconnection site migrates northward at graduahy 
lower speeds, longer and slightly less sheared loops are 
formed one above the other, and this results in the in- 
ferred increase of the altitude of the LT source and the 
separation of the FPs from the NL. 

Phase II is in agreement with the classical CHSKP 
picture of two-ribbon flares, while Phase I is not. One 
possible physical explanation for Phase I^ suggested by 
iJi et alJ (2007) is the magnetic "implosion" conjecture 
(|Hudsonl 2000) that predicts contraction of field lines 
durin g a fiare as a c onsequence of explosive energy re- 
lease. [11^1^ (|2007l ) further found that a magnetic field 
with shorter, lower lying, and less sheared field lines in- 
deed contains less free energy. N ote that, in or der to 
explain the LT descending motion, IVeronig et al.l (|2006) 
proposed a collapsing magnetic trap model, which, how- 
ever, cannot explain the initial FP motion toward one 
another para llel to the NL. In contrast, the scenario of 
ISomov et al.l (2002) for the approaching FPs does not 
predict the decrease of the LT height. 

4. TEMPORAL CORRELATIONS OF CONJUGATE 
FOOTPOINTS 

We now examine the temporal evolution of and corre- 
lations between various quantities of the two conjugate 
FPs, particularly spectral, spatial, and magnetic field pa- 
rameters, which are summarized in Table [2l 

4.1. Spectral Correlations 

Figure [9]a shows the history of the photon fiuxes of 
E-FP (/i, blue diamonds) and W-FP (/2, red crosses) 
at 50 keV obtained from the power-law fits in the 50- 
150 keV range mentioned in § 12. 2[ We find that the two 
fiuxes follow each other closely in their temporal trends 
and E-FP is always brighter than W-FP except for the 
first time interval. The correlation of the fiuxes can also 
be seen in Figure [TQl a where one fiux is plotted vs. the 
other. A linear regression is shown as the thick dashed 
line and given in Table [2l The correlation coefficients 
listed in Table [2] indicate a very high correlation in ei- 
ther a linear or a nonlinear sense. Such a correlation 
is expected for conjugate HXR FPs, since they are be- 
lieved to be produced by similar populations of nonther- 
mal electrons that escape the same accelerati on region 
(belie ved to be at/near the LT source; Petrosi an fc Liul 
l2004t iLiu et al.l 120081) and travel down opposite legs of 
the same magnetic loop to reach the chromosphere. 

We show the corresponding power-law indexes (7) of 
the two FPs vs. time in Figure [9]ri and one index vs. the 
other in Figure \Wb . Again we find that the two indexes 
are closely correlated, as can be seen from the large cor- 
relation coefficients (Table [2j). The E-FP spectrum, how- 
ever, is slightly harder than the W-FP spectrum, which 
it is persistent most of the time. The results from long in- 
tegration intervals (2-3 minutes, not shown), which have 
better count statistics, exhibit the same pattern. We av- 
eraged the index values of the six 2 minute intervals cov- 
ering 20:40:40-20:52:40 UT, after which the uncertainties 

^ Considering the presence of multiple FPs (see Fig. [3]), this 
phase might prov ide evidence for the tether-cutting model of 
IMoore et all (|200l[ ). 



become large due to low count rates. This average gives 
(71) = 3.63 ± 0.06 for E-FP and (72) = 3.79 ± 0.11 for 
W-FP. Their difference of (72) - (71) = 0.15 ± 0.13 is 
marginally significant at the la level. 

Let us compare the HXR fluxes and spectral indexes 
of the two FPs. As can be seen in Figures [9]a and [9]ri, 
during HXR Peak 1 (before ^2=20:48 UT), the fluxes 
and indexes are anti-correlated, i.e., they follow the gen- 
eral ''soft-hard-soft " (SHS) tren d obser ved in many other 
flares (e.g., iParks fc Winckle^ [19691 : IKane fc Ande rsonl 
I197QI)- However, during Peak 2 (after ^2), the indexes 
decrease through the HXR maximum and then vary only 
slightly (with relatively large uncertainties) around a 
constant level of 3.0. This trend can be characterized 
as "soft-hard-hard" (SHH). This flux-index relationship 
can also be seen in Figure [TOfc where the index averaged 
between the two FPs is plotted against the average flux. 
Note that the spectral index values during the late de- 
clining phase of the flare are even smaller than those at 
the maximum of the main HXR Peak 1. In this sense, the 
overall spectral variation can be characterized as "soft- 
hard-soft-harder" . As we noted earlier, there were ener- 
getic protons detected in interplanetary space by GOES 
and ACE following the flare. These observations, when 
taken together, are consistent with the statistical result 
of Kiphnger (1995) that this type of flare with progres- 
sive spectral hardening tends to be a ssociated with SEP 
events (also see ^Saldanha et al.ll2008l ). As we also noted, 
stron g ^amma-ray line em ission was detected during this 
flare (Hurford et al."2006) , which indicates a signiflcantly 
large population of accelerated protons at the Sun, but 
the relation to the SEPs at 1 AU is unclear. 

4.2. Spatial Correlations 

We now examine the spatial evolution of two FPs. In 
§ [3] we focused on their relative motion, while here we 
compare their individual motions. Figure [8]e shows the 
perpendicular distance of each FP from the north-south 
NL {red dashed^ Fig. [Th) as a function of time. Linear 
flts of the full flare duration indicate mean velocities of 
{vi_) = 36 ± 1 for E-FP and 13 ± 1km s'^ f or W-FP 
These velocities are similar to those found bv IXu et al.l 
(|2004l ) for near infrared ribbons in this event. We also 
calculated the total velocities of the FP centroids, i.e., 
V = {v'^ + 'Uy )^/^, where v± and are the components 
perpendicular and parallel to the NL, respectively. The 
two resulting velocities have a linear temporal correla- 
tion at a 5<j level (see Table [2j), which again provides 
evidence of the causal connection between the conjugate 
FPs. However, the individual component, v_\_ or , alone 
does not exhibit any noticeable correlation between the 
two FPs. 

Figure [8]a shows the distances from the LT centroid 
to the centroids of E-FP (/i, diamonds) and W-FP (^2, 
crosses) along the model semi-circular loop (see, e.g.. 
Fig. [3)) as a function of time. Each curve follows the 
same general increase as the corresponding distance from 
the NL shown in Figure [8]e, but the distance to E-FP is 
smaller than that to W-FP most of the time. We esti- 
mated the coronal column densities from the LT source 
to the transition region at the two FPs (see Appendix [Pl 
for details) to be A/'tr,i = 0.5nLT(^i — ^), where z=l, 2, 
using the distances Z^, LT density nur and equivalent ra- 
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TABLE 2 

Correlation coefficients and linear regressions between various parameters of the conjugate footpoints for the full 
FLARE DURATION (20:40:40-20:59:40 UT). Magnetic field strengths are in units of 100 Gauss. 





Subscripts: 




#of 




signif. 


Linear regression (between 


1: 


E-FP, 2: W-FP 




cr's 






quantities in first two columns) 


h 


h 


0.98 


8 


0.97 


10-^^ 


/2 = (-0.3 ± 0.1) + (0.41 =b 0.01)/i 


71 


72 


0.90 


7 


0.89 




72 = (-0.5 =b 0.3) + (1.17 zb 0.07)71 


Bi 


B2 


0.39 


3 


0.40 


10"^ 


B2 = (-9.2 ± 0.8) + (3.7 ± 0.1)5i 


Bi 


logio ^1 


0.50 


4 


0.49 


10-^ 


h = (0.012 ± 0.003) X io(o-55±o.o2)Si 


B2 


loglO ^2 


0.82 


6 


0.84 


10-^^ 


/2 = (0.19 ± 0.01) xl0(°-^3"=^°-°°2)^2 


B 


logio ^ 


0.77 


6 


0.84 




(0.079 ±0.007) X io(o-252±o.oo5)s 


Vl 


V2 


0.63 


5 


0.29 


10"^ 


V2 = (-93± 10) + (2.8±0.1)^;i 



Note. — r^: Pearson's linear correlation coefficient; 
# of cr's: multiple of la uncertainty of r^, where a = 1/(57)^^^ 
intervals) ; 

rs and signif.: Spearman's rank correlation coefficient and significance level; 
/ (photons s-^ cm-2 keV"^): HXR flux at 50 keV; 

7: spectral index between 50 and 150 keV (regression done only for 20:40:40-20:52:40 UT); 
B (in 100 G): magnetic field; 

B (in 100 G) and / (photons s~^ cm~^ keV~^): B and / averaged between E-FP and W-FP 



0.13, in which 57 is the number of data points (or time 



V (km 



FP velocity. 



(0): FP HXR flux I ot 50 keV 
■ E-FP o 
2: W-FP X - 




(e)| coronal column density 
to transition region at FP- 



c^j 20:44 20:48 20:52 20:56 21:00 

Time (UT, 2003/10/29) 

Fig. 9. — Time profiles of X-ray and magnetic field parameters, 
(a) HXR fluxes at 50 keV of E-FP {blue diamonds) and W-FP {red 
crosses) obtained from power-law fits in the 50-150 keV range. 
The same color and symbol convention for E- and W-FP holds 
for the other panels. The vertical dotted lines correspond to the 
transitional times of ti and t2 as shown in Fig. [2] and ts, the 
maximum of HXR Peak 2. {b) SOHO MDI magnetic field strengths 
registered at the two FPs. (c) Ratios of the 50 keV fluxes (W-to-E) 
and magnetic fields (E-to-W) of the two FPs. {d) HXR spectral 
indexes of the two FPs from the same fits as in (a). The arrow 
marks the end of the s ix 2 minute integration intervals for averaging 
the index values (see ^4.1|) . (e) Estimated coronal column densities 
from the edge of the LT source to the transition region at the two 
FPs. 
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Fig. 10. — Various correlations between the HXR and magnetic 
field parameters of the two FPs. (a) HXR flux at 50 keV (Fig.[9h) of 
W-FP vs. E-FP. The dashed line here and in other panels represents 
a linear regression of the data. The color bar indicates the common 
time evolution for all panels of this figure. (6) Same as (a) but for 
the spectral indexes (7) shown in Fig. [9ji. The linear regression is 
done for the data points up to 20:52:40 UT. (c) Average (between 
the two FPs) spectral index 7 vs. average flux / at 50 keV. The 
two open circles here and in {d) mark the times t2 and ts shown in 
Fig. [9] A "soft-hard-soft" (SHS) variation is present early in the 
flare, {d) Same as (a) but for the magnetic field strengths of the 
two FPs as shown in Fig. [9^. (e) Logarithmic HXR flux at 50 keV 
vs. the corresponding magnetic field strength for E-FP {diamonds) 
and W-FP {plus signs), {j) Same as (e) but for the average values. 



Asymmetric Conjugate HXR Footpoints 



11 



dius r obtained earlier (Fig. [6]c). The results in Figure [9]e 
show that there is a large relative difference from 20:44 
to 20:48 UT during HXR Peak 1 but a smaller difference 
during Peak 2. Implications of these different column 
densities will be addressed in § 15.21 and § 15.31 

4.3. Magnetic Field Correlation 

The magnetic field strengths of the two FPs were ob- 
tained from SOHO MDI magnetograms (e.g., see Fig.[7la) 
through the following steps: (1) We first selected a pre- 
flare^ magnetogram at 08:30:35 UT and coaligned it with 
the RHESSI pointing and field of view using the offsets 
found in Appendix[Bl (2) For each time interval, the 90% 
brightness contour (not necessarily a resolved source) of 
each RHESSI FP source was rotated back to the corre- 
sponding position at the time of the MDI map to account 
for the solar rotation. Then the MDI pixels enclosed in 
this contour were averaged to give a value of the mag- 
netic field for this FP, and the standard deviation of these 
pixels combined with the nominal 20 G MDI noise was 
used as the uncertainty. (3) The above two steps were re- 
peated for each of the ten MDI magnetograms recorded 
between 20:25 and 20:35 UT at a one minute cadence 
(excluding the one at 20:28 UT that are contaminated 
by artificial pixel spikes). The average of the ten inde- 
pendent measurements was used as the final result for 
the magnetic field (Fig. [9]?)), and error propagation gave 
the final uncertainty (~5-10%). 

As shown in Figure [9]?), the magnetic field strength of 
W-FP generally decreases with time, while that of E-FP 
fluctuates about its mean value. Most of the time (es- 
pecially before ^3=20:51 UT), the W-FP field is stronger 
than the E-FP field, while their fractional difference gen- 
erally decreases as time proceeds. The temporal varia- 
tions of the two field strengths are only weakly corre- 
lated (again, particularly before t^) at the 3cr level (see 
Table [2j), as can also be seen in Figure [TOb. 

4.4. Inter- correlations Among Spectral, Spatial, and 
Magnetic Field Parameters 

Here we check the relationship between the HXR fluxes 
and the magnetic fields of the two conjugate FPs. We 
plot the logarithmic HXR flux vs. the magnetic field 
strength for each FP in Figure [TOl e. As we can see, 
the flux is correlated with the field strength for each 
FP (see Table [2] for the correlation coefficients and lin- 
ear regressions). The logarithmic average flux (/) and 
magnetic field {B) of the two FPs are shown one vs. 
the other in Figure [TOjf. A linear relationship, as shown 
by the thick dashed line with a correlation coefficient of 
= 0.77 ± 0.13, is clearly present. In other words, / 
is exponentially {nonlinearly) correlated with 5, the ex- 
pression of which is listed in Table O Since / is correlated 
with 5, the "soft-hard-soft" type of relationship between 
7 and / during the early phase (before ^3=20:51 UT) 

^ Note that magnetograms during the flare cannot be used due to 
temporary artiflcial changes in the measured fleld strength (see Ap- 
pendix EI- Permanent real changes have been observed before and 
after many X-class flares (e.g., Wang et al. 2002; Sudol &; Harvev 
|2005), and have been interpreted as magnetic fleld changes in direc- 
tion rather than in strength (Sudol & Harvev 2005). Consequently 
we use the preflare fleld as the best approximation available to us 
for the fleld during the flare. 



translates to that between 7 and B. Namely, 7 is anti- 
correlated with B. 

Finally we check the relationship of the apparent mo- 
tions and magnetic fields of the two FPs. As noted above, 
E-FP moves faster than W-FP away from (perpendic- 
ular to) the magnetic NL, while E-FP is located in a 
weaker magnetic field. This relationship means that, as 
expected, about the same amount of magnetic flux is an- 
nihilated from each polarity, since B{v±) is proportional 
to the magnetic reconnection rate. However, the mag- 
netic fluxes swept by the two FPs, / B{v±)dt integrated 
over the full flare duration, differ b y 44% of their aver- 
age va lue. This is not surprising, as lFletcher fc HudsonI 
(|200lf ) found a similar flux mismatch and offered various 
explanations. 

4.5. Discussion on Implications of Various Correlations 

The above correlation (Fig. [TOjf ) between the average 
HXR flux and magnetic field strength reveals important 
information about the magnetic reconnection and par- 
ticle acceleration processes. Here we speculate on two 
alternative possible scenarios in terms of contemporary 
acceleration theories. 

1. The nonlinear (exponential) nature of the I-B cor- 
relation suggests that particle acceleration is very sensi- 
tive to the magnetic field strength, if we assume that the 
measured photospheric field strengths scale with that in 
the coronal a cceleration region. The stochastic accelera- 
tion model of lPetrosian fc Liul (|2004l ) offers the following 
predictions: (1) The level of turbulence that determines 
the number of accelerated electrons is proportional to 
SB'^^ where SB is the magnetic field amplitude of plasma 
waves. The acceleration rate that determines the spec- 
tral hardness of accelerated electrons is proportional to 
SB^/B. An increasing field strength B will result in an 
increasing flux and spectral hardness of accelerated elec- 
trons, if both 5B'^ and 6B^/B also increase. (2) The 
relative efficiency of acceleration of electrons and thus 
their spectral hardness increase with decreasing values of 
the ratio of electron plasma frequency to gyro-frequency, 
a = uOpe/Vte (X 1/B. These predictions are qualitatively 
consistent with the observations that the magnetic field 
strength correlates with the HXR flux and anti-correlates 
with the spectral index. 

2. Alternatively, noting the roughly constant veloci- 
ties {v±^ Fig. [8]e) of the two FPs perpendicular to the 
magnetic NL, the above I-B correlation simply trans- 
lates into the correlation between the HXR production 
rate and the magnetic flux annihilation rate or reconnec- 
tion rate, Bv±. Furthermore, since Bv± is believed to 
be proportional to the electric field in the reconnection 
region (Forbes fc Linll2000l : iQiu et al.ll2002h . it then fol- 
lows that the particle acceleration rate correlates with 
the electric field. Accord i ng; to the e lectric field accel- 
eration model of iHolmanI (|l985l ) and i Benka fc HolmanI 
(1994), a larger electric field results in a larger high- 
energy cutoff (^max) for the electr on spectrum, which 
can lead to a harder HXR spectrum (|Holmanll2003l ). This 
is consistent with the observed anti-correlation between 
the magnetic field strength and spectral index. Note that 
in the classical (Petschek 1964) model, the small cross- 
section of the current sheet cannot account for the typi- 
cally large f lux of accelerated electrons of 10^^-10^^ elec- 
trons s~^ (iMiller et al.lll997l ). i.e., the so-called "number 
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problem". However, obs ervations and mag net ohydrody- 
namic simulations (e.g. jKliem et al.ll2QQCl ) have implied 
that the reconnecting current sheet involves small-scale 
electric fields around multiple, spatially separated mag- 
netic X- a nd /or -points in a fragmented topology (see a 
review bv I Aschwanden 2002) . We speculate that, when 
particle acceleration takes place in such a fragmented 
current sheet, the numbe r problem can be amelior ated 
(also see a discussion bv iHannah fc Fletcherl [2"006f ) . In 
this case, the above discussion remains valid if the small- 
scale electric fields are scaled with the macroscopic po- 
tential dr op and thus with B v±. For comparison, we 
note that iKrucker et al.l (|2005l ) studied the motion of E- 
FP alone in this fiare and also found a rough temporal 
correlation between the HXR fiux and reconnection rate, 
represented by Bv or B^v, where v includes the velocity 
both perpendicular and parallel to the NL. 

5. HARD X-RAY FOOTPOINT ASYMMETRIES 



As mentioned above and partly noted by jLiu et al. 
(j2004c[ ). IXu et all (j2004l ). and Krucker^^^ pOOS), the 
two conjugate FPs exhibit the following asymmetric 
characteristics: (1) the brighter E-FP is located in a 
weaker^ negative magnetic field, while the dimmer W- 
FP is located in a stronger, positive field; (2) the two 
FPs have very similar spectral shapes with E-FP being 
slightly harder; (3) E-FP is located closer to the LT than 
W-FP; (4) E-FP moves faster away from the magnetic 
NL than W-F. These asymmetries are summarized in 
Tables 

We explore in this section different possibilities that 
can cause such asymmetries, particularly the asymmet- 
ric HXR fiuxes and spectra. Various physical processes 
can contribute and they fall into two categories ac- 
cording to their origins: (1) asymmetry during parti- 
cle acceleration, and (2) asymmetry arising from parti- 
cle transport. The second category includes effects of 
magnetic mirroring and column density, which will be 
examined in what follows ( §5.11[5^ . Other transport 
e ffects and the first category will be discussed later in 
§ 15.41 We use both th e fiux ratio Rr = h/h and 
the asymm etry (c.f., Asch wanden et al] 1 19991 ) defined by 
[Alexander fc Metcalf (2002h : 

A = {h- h)/{h + h) = (1 - Ri)/(l + Ri). (1) 

to quantify the asymmetric HXR fiuxes, with A = ±1 be- 
ing 100% asymmetry and A = being perfect symmetry. 



5.1. Magnetic Mirroring 

Asymmetric magnetic mirroring is commonly cited to 
explain asymmetric HXR fiuxes observed at conjugate 
FPs. We examine to what extent mirroring alone can 
explain the observations of this fiare. For simplicity, we 
make the following assumptions for our analysis below: 

(1) Disregard all non-adiabatic effects of particle trans- 
port, i.e., energy losses and pitch-angle diffusion due to 
Coulomb collisions. By this assumption, the magnetic 
moment of a particle is conserved and mirroring is the 
only effect that changes the pitch angle when the particle 
travels in the loop and outside the acceleration region. 

(2) Assume an isotropic pitch-angle distribution of the 
electrons at all energies upon release from the accelera- 
tion region. (3) Disregard details of bremsstrahlung, and 



TABLE 3 

Asymmetric characteristics of the conjugate footpoints 
(e-fp and w-fp): mean, median, and their e-to-w ratio of 
various quantities. 



Mean 



Median 





E 


W E/W 


E 


W E/W 


/ 


13.7 


6.1 2.2 


8.9 


4.6 1.9 


B 


520 


960 0.55 


520 


1010 0.51 


7 


3.63 


3.79 0.96 


3.4 


3.5 0.96 


/ 


41 


60 0.69 


38 


60 0.64 


iVtr 


1.2 


2.1 0.60 


1.2 


2.5 0.50 




36 


13 2.8 







Note. — / (photons cm~^ keV"^): HXR flux at 
50 keV; 

B (Gauss): magnetic field strength; 

7: spectral index between 50 and 150 keV; 

/ (arcsecs): distance from the LT centroid to the FP centroid 

along the semi-circular loop; 

A^tr (10^° cm~^): coronal column density from the edge of the 
LT source to the transition region at the FP. 
v± (km s~^) : velocity perpendicular to the NL; 
The mean spectral indexes are the averages of the results of 
the six 2 minute integration intervals between 20:40:40 and 
20:52:40 UT (see ^iTTI) . while all the other values listed here 
are from the results of the 57 short intervals throughout the 
flare. 

assume that the nonthermal HXR flux is proportional to 
the precipitating electron flux at the FP.'^ 

The loss-cone angle for magnetic mirroring is given 
as Oi = arcsin(5o/5^)^/^ (i=l for E-FP, 2 for W-FP), 
where Bq is the magnetic field strength at the injection 
site in the corona where particles escape from the ac- 
celeration region, and Bi is the field strength at the 
ith FP in the chromosphere. By the isotropy assump- 
tion, the fractional flux of the forward moving elec- 
trons that will directly precipitate to the chromosphere 
(whose pitch angle is located inside the loss cone) can 
be evaluated by integrating over the solid angle (also see 
[Alexander fc Metcalf 2002 ): 



Fi = — [ dn = — [ dcj) j \mede = 1- in, 

27r J 27r Jq Jq 



(2) 



where the pitch- angle cosine is jUi = cos^i = (1 — 
Bo/Bi)^^'^. If there is strong mirroring, i.e., Bq <^ Bi, 
we have /i^ 1 — Bo/{2Bi), and if there is no mirroring, 
i.e., ^0 = Bi, we have /i^ = 1. By our assumptions (1) 
and (3) above, such a fraction should be independent of 
electron energy and is proportional to the HXR flux li 
at the corresponding FP. It then follows that 




ifBo^Bi,{i = l, 2), 

if Bq = Bi, Bq <^ B2, 



(3) 



(4) 



the second case of which corresponds to the possibility 
that mirroring occurs only at one FP, but the required 
condition Bi <C B2 does not apply to this flare. In either 

^ This flux includes the precipitation of electrons previously re- 
flected by mirroring back to the acceleration region at the LT where 
they may be scattered and/or re- accelerated, presumably by tur- 
bulence. 
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case, the HXR flux ratio should be correlated with the 
inverse of the magnetic field ratio, R^^. This result is 
co nsistent with that o f the strong diffusion case obtained 
bv lMelrose fc White! (|l98ir ). 

As shown in Table [31 the mean/median HXR fiux of E- 
FP is about twice that of W-FP, while the mean/median 
magnetic field strength of E-FP is about a factor of two 
smaller. This is consistent with the mirroring effect in 
the average sense of the whole flare duration. We can 
check if this relationship also holds at different times. 
Figure [9]c shows the W-to-E ratio {Rj = h/h) of the 
HXR fluxes and the E-to-W ratio (i?^^ = B1/B2) of the 
fleld strengths of the two FPs as a function of time. We 
flnd a temporal correlation between the two ratios during 
the flrst ~3 minutes when both flrst decrease and then 
increase. In the middle stage (20:43-20:51 UT) of the 
flare, both ratios remain roughly constant with marginal 
fluctuations and similar mean values of {Ri) = 0.43 and 
{R-^^) = 0.50. After 20:51 UT, the magnetic fleld ratio 
increases signiflcantly with large fluctuations and exceeds 
unity in 6 of the 57 time intervals, while the flux ratio 
remains at about the same level as before. The behav- 
ior of the two ratios before 20:51 UT is expected from 
magnetic mirroring, but their signiflcant difference after 
20:51 UT cannot be explained by mirroring alone. 

5.2. Column Density 

Another transport effect that can cause asymmetric 
HXR FPs is different coronal column densities experi- 
enced by electrons in traveling from the acceleration re- 
gion t o the tra nsition region at the two FPs (^Emslie et al.l 
120031 : [Liul 12006. . p. 73). The effective column density is 
^tr,eff = ^tr/ (m) , whcrc (/i) is the average pitch angle co- 
sine, and A^tr = j^^^ n{s)ds is the coronal column density 
to the transition region at distance s = 5tr along the mag- 
netic fleld line with n{s) being the ambient electron num- 
ber density. A difference in (/i), 5tr, and/or n{s) between 
the two legs of the flare loop can lead to different effective 
column densities. (1) Different pitch- angle distributions 
can be caused by asymmetric magnetic mirroring and/or 
asymmetric acceleration. (2) Different path lengths 5tr 
can be caused by a magnetic re connection site located 
away from the middle of the loop (|Falewicz fc Siarkowskil 
[2007) . (3) Different densities n{s) can also occur be- 
cause magnetic reconnection takes place between fleld 
lines that are previously not connected and their associ- 
ated densities are not necessarily the same. It takes time 
(on the order of the sound travel time, > tens of sec- 
onds) for the newly reconnected loop to reach a density 
equilibrium, but the observed HXRs could be produced 
before then. 

Column density asymmetry affects the FP asymmetry 
in two ways, since both energy losses and pitch-angle 
scattering due to Coulomb collisions take place at about 
the same rate that is proportional to the column density: 
(1) Column density asymmetry is related to energy losses 
and the way we calculate the FP photon flux in § 12.21 (in- 
tegrating HXR photons primarily produced below the 
transition region). Electrons with an initial energy of 
E are stopped after traveling through a column density 
A^stop(cm-2)/ ^ 10i^[^(keV)]2. If A^^top is smaller than 
or comparable to the column density to the transition 
region (A/tr), in one half of the loop with a larger column 



density, there are more electrons stopped in the leg and 
thus less electrons reaching the transition region. This 
results in more HXRs produced in the leg and less HXRs 
beneath the transition region (counted as the FP flux). 
(2) Different Coulomb scattering rates result from differ- 
ent column densities on the two sides of the loop, which 
can cause different pitch- angle distributions, even if the 
particles are injected with symmetrical pitch angles from 
the acceleration region. 

Focusing on the energy dependence of FP HXR asym- 
metry, we present below an estimate of the column den- 
sity effect alone^ while assuming no magnetic mirroring 
and identical pitch- angle distributions (same (/i)) in the 
two loop legs. The relevant formulae are derived in Ap- 
pendix O We further assumed that identical power-law 
electron fluxes with a spectral index 5 are injected into 
the two legs of the flare loop, which have the same am- 
bient density but different path lengths 5tr to the FPs. 
Our goal is to examine if this scenario can yield photon 
fluxes and spectra consistent with observations for both 
FPs. For each of the 57 time intervals shown in Figure [9l 
we flrst used the E-FP column density Ntr (see §4.21 and 
Fig. [9]e) from the edge of the LT source (assumed to be 
the acceleration region) to obtain its dimensionless form 
Ttr = Ntr/{^ X lO^^cm"^), which ranges from 5 x 10~^ 
to 1 X 10~^. We then substituted rtr into 

/ L , 1 \ 1-5/2 

/Fp(rtr, k) = Aok-^'-'^ ^1 + nr^^j , (5) 

rewritten from equation ()C4p . where Aq is the normaliza- 
tion for the thick-target flux and k is the photon energy 
in units of the rest electron energy 511 keV. With this 
equation, we fltted the E-FP spectrum above 50 keV in 
a least-squares sense by iteratively adjusting the free pa- 
rameters Aq and S. Using the resulting Aq and 5 and 
W-FP's Ttr, we then calculated the W-FP spectru m b y 
equation ([5]) and the W-to-E flux ratio by equation ()C5p . 

Figure [TTb shows an example of the spectra of the two 
FPs and their model predictions, together with the cor- 
responding thick-target spectrum produced by the same 
power-law electron flux. We only trust the observed FP 
spectra >50keV due to pileup, as noted earlier. As ex- 
pected, the model FP fluxes are reduced from the thick- 
target flux, especially at low energies, because low-energy 
electrons are more susceptible to collisional energy loss 
and pitch-angle scattering. This results in a spectral 
flattening (hardening) in the FP X-ray spectrum. Be- 
cause of its larger column density, the W-FP's model 
spectrum exhibits more flux reduction at a given energy 
and a flattening to a higher energy. Above 50 keV the 
model spectrum of the brighter E-FP flts the data very 
well. However, that of the dimmer W-FP does not flt the 
data at all, since the model flux is much greater (e.g., at 
77.5 keV / = 4.0 vs. 2.6 photons s"^ cm"^ keV"^) and 
harder (7 = 3.5 vs. 4.1) than the observed flux, and even 
harder than the E-FP flux (7 = 3.5 vs. 3.9). This can 
be best seen in Figure [TTb that shows the data {plus 
signs) and model {solid line) ratios of the W-to-E FP 
flux. The data ratio generally decreases with energy or 
stays roughly constant above 50keV within uncertainties, 
but the model ratio is an increasing function of energy. 
These trends generally hold throughout the flare as can 
be seen from the history of HXR fluxes and spectral in- 
dexes shown in Figures [TTb and[TT1d. 
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Fig. 11. — Effects of asymmetric coronal column densities, (a) Photon fluxes of E-FP (diamonds) and W-FP (crosses) vs. energy at 
20:44:40-20:45:00 UT (same as Fig. [Sja), superimposed with model fluxes (lines) evaluated l3elow the transition region using eq. (|5]). The 
solid line shows the power-law thick-target flux with an index of 7thick =4.1. The legend shows the fluxes (/) at a reference energy of 
77.5 keV and the spectral indexes (7) above 50 keV of the two FPs for both the data a nd m odel, (b) Ratio of the W-to-E FP fluxes shown 
in (a), with the plus signs for the data and the solid line for the model (given by eq. IC5)) . (c) Observed (symbols) and modeled (lines) 
fluxes of E- and W-FP at 77.5 keV vs. time (cf. Fig. [9la at 50 keV). (d) Same as (c) but for the spectral indexes above 50 keV. The data is 
the same as in Fig. ^ but without error bars for clarity. 



In summary, the model predicts a much harder photon 
spectrum for the dimmer W-FP with the larger coronal 
column density, while according to the observations the 
dimmer W-FP is as hard as or slightly s ofter than the 
brighter E-FP (see Fig.^). Saint-Hilaire et al. (2008) 
reported similar results that the majority of the brighter 
FPs in their 172 pairs of FPs during the HXR peaks of 
53 flares tend to have harder spectra. In addition, the 
differences between the model HXR fluxes of the two FPs 
are too small to explain the observations (Fig.fTTfc). One 
may attempt to increase the difference between the col- 
umn densities in order to increase the flux difference and 
thus to merge this gap between the model and data, but 
the discrepancy of the spectral indexes would be exacer- 
bated since the W-FP spectrum would be relatively even 
harder. Therefore, we conclude that the column density 
effect alone cannot provide a self-consistent explanation 
for both the HXR fluxes and spectra of the asymmetric 
F Ps observed here. 

iFalewicz fc Siarkowskil (|20Q7[ ). however, found in three 
flares that the HXR flux ratios of asymmetric FPs were 
consistent (within a factor of 2) with the predictions from 
asymmetric column densities. While that scenario may 
apply to those flares, we should note that their broad 
band Yohkoh observations set less stringent constraints 
than our high resolution RHESSI observations, which 
can lead to different conclusions. In particular, their 
analyses were limited to images in the Ml (23-33 keV) 
and M2 (33-55 keV) bands, where the column density 
effect is more pronounced and contaminated from ther- 
mal emission is possible, while our observations cover a 
higher and wider energy range of 50-150 keV. 



5.3. Magnetic Mirroring and Column Density Combined 



We have seen from the above discussions that each of 
the two transport effects alone can only explain to some 
extent the observed FP asymmetries: (1) Asymmetric 
magnetic mirroring is consistent with the asymmetric 
HXR fluxes in the average sense of the flare duration, 
but it has difliculties in accounting for the flux asym- 
metry later in the flare. (2) Asymmetric column den- 
sities in the two legs of the flare loop are qualitatively 
consistent with the asymmetric HXR fluxes, but their 
quantitative predictions of fluxes and spectral hardness 
contradict the observations. These two transport effects, 
in reality, operate at the same time^ because electrons 
experience Coulomb collisions while being mirrored back 
and forth in the loop, and thus the collisionless (adia- 
batic) assumption that we adopted earlier for simplicity 
for magnetic mirroring needs to be dropped. In particu- 
lar, since W-FP has stronger mirroring (than E-FP), the 
average pitch angle of electrons impinging there is larger, 
and thus the effective column density A/tr,eff = Mr/(/^) 
is greater than previously thought. This can enhance 
the column density asymmetry. In what follows, we at- 
tempt to provide an explanation for some aspects of the 
observations by combining the two effects. 

From the above discussion and the observations pre- 
sented in § m we should pay attention to the distinction 
between the two HXR peaks. As shown in Figure [9]c, 
during Peak 1 (< ^2=20:48 UT) the FP HXR flux asym- 
metry seems to be mainly controlled by magnetic mirror- 
ing, while during Peak 2 (> ^2), especially after the HXR 
maximum at ^3=20:51 UT, this control seemingly fails. 
A viable explanation for the two-peak distinction is that 
(1) at early times, the densities (Fig. [6^) and lengths 
(Fig. [8]a) of the loops are small, resulting in small coro- 
nal column densities (Fig. [9]e) from the acceleration site 
to the FPs. Energy losses and pitch-angle scattering due 
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to Coulomb collisions are less important, and therefore 
the rates of electron precipitation to the FPs are mainly 
governed by mirroring. (2) Later in the flare, as the 
loop densities and lengths have increased considerably, 
the column densities become larger and the collisional 
effects become more important than before in shaping 
the observed FP flux asymmetry. In addition, since mag- 
netic mirroring depends on the gradient dlnB/dN (e.g., 
iLeach fc Petrosianlll98ll ). this effect becomes less impor- 
tant when the column density N increases faster than the 
relative change of magnetic field from the LT to the FP, 
which is possibly the case later in the flare. Therefore, 
at later times, the prediction of magnetic mirroring alone 
tends to deviate from the data, which might be explained 
by the two transport effects combined. The outcome of 
the combination can be modeled quantitatively as dis- 
cussed later in § [6] for future work. 

There are several coincidences with the two-peak di- 
vision which seem to have causal connections: (1) As 
noted in § 14. If the correlation between the HXR fluxes 
and spectral indexes (Figs. [9]a and[9H) can be described 
as common "soft-hard-soft" during Peak 1 and as "soft- 
hard-hard" during Peak 2. The spectral hardening at 
later times may be associated with the increasing loop 
column densities (Fig. [9]e), due to collision-caused hard- 
ening mentioned above ( §5.2p . (2) During Peak 1, the 
magnetic fields at the two FPs are weakly correlated with 
each other (Fig. [9?)), while this correlation becomes pro- 
gressively nonexistent during Peak 2, especially after its 
maximum ^3=20:51 UT, possibly because of longer loops. 
(3) The transition (^2 =20:48 UT) between the two HXR 
peaks coincides with the sudden jump in the positions of 
both FPs (Figs. [Th and[8]e), the dip in the loop length 
(Fig. [8]a), and the valley in the magnetic field strengths 
(Fig. 13)). This points to the start of the new episode of 
energy release of Peak 2, presumably associated with a 
new series of loops that have physical conditions different 
from those in Peak 1. This transition may be related to 
the different behaviors of magnetic mirroring during the 
two peaks noted above. 

5.4. Discussion on Other A symmetry- causing Effects 

Here we briefly discuss transport effects other than 
magnetic mirroring and column density, and acceleration 
related effects that can contribute to the observed HXR 
flux and spectral asymmetries. 

Asymmetry during the transport process — (1) Non- 
uniform target ionization: In the above analy- 
sis we assumed a fully ionized background in the 
path of high-energy electrons, but in reality the back- 
ground targets vary from fully ionized in the corona 
to neutral in the chromosphere. The presence of 
neutral atoms reduces the rates of long-range colli- 
sional energy losses and pitch-angle scatte ring; and thus 
increases the bremsstra h lun^ efficiency (iBrownl Il973l : 
Heach & PetrosianI Il98ll : iKontar et aP 120021 ). In this 
flare, at the E-FP with weaker magnetic mirroring elec- 
trons can penetrate deeper into the chromosphere and 
thus encounter more neutral atoms. This may result in a 
higher HXR flux and harder spectrum in the 50-150 keV 
range at E-FP than at W-FP, qualitatively consistent 
with the observations. (2) Relativistic beaming and 
photospheric albedo: At the E-FP with weaker mir- 
roring, the angular distribution of electrons are more con- 



centrated to the forward direction down to the photo- 
sphere. When the FPs are seen on the solar disk from 
above, the increasing importanc e with energy of the for- 
ward relativistic beaming effect (i McTiernan fc PetrosianI 
1991) results in a smaller fraction of high-energy pho- 
tons emitted upward at E-FP. Meanwhile, since relatively 
more photons are (beamed) emitted downward at E-FP, 
albedo or Compton back-s cattering (|Langer fc PetrosianI 
119771 : iBai fc RamatvlfT978h is stronger there. Both effects 
can cause a softer X-ray spectrum at E-FP than at W- 
FP, competing with other effects mentioned above, which 
may explain why the spectral indexes are so close ( (72 ) — 
(71) = 0.15 ±0.13). (3) Return currents and the asso- 
ciated electric field decrease the energy of the downward- 
streaming electrons, with the major impact being on th e 
lower-energy electrons (Zharkova fc Gordovskvvl l2006l ). 
Different precipitating electron beam fluxes in the two 
legs of the flare loop may induce different return current 
densities, and thus can result in different HXR fluxes and 
spectral shapes at the two FPs. 

Intrinsic asymmetry arising from the particle accel- 
eration process — An energy-dependent FP HXR flux 
asymmetry (|A|, eq. [l]), which has a maximum in 
the intermediate energy range (20-40 keV) and de- 
creases toward both low and high energies, was found 
by Alexander fc Metcah (2002). This was attributed 
to an asymmetric, energy-dependent accelerator from 
which more electrons are inje cted preferentially into one 
of th e two legs of the loop (|McClements fc Alexanderl 
I2OO5I ). Asymmetric electron beams can be produced 
by the electric field in a rec onnection current sheet 
([Zharkova fc Gordovskvvl l2004l ). For the flare under 
study no reliable asymmetry can be obtained below 
50 keV due to pulse pileup; above 50 keV asymmetry ei- 
ther increases with energy (see the flux ratio in Fig. [5]??) 
or remains constant, opposite to the decreasing asymme- 
try reported by Alexan der fc Metcalfl_( [2002) in this en- 
ergy r ange. From Yohkoh HXT data ^Aschw anden et al.l 
(|l999l ) also found no general energy-dependent pattern of 
flux asymmetry. Furthermore, if acceleration by plasma 
waves is the dominant mechanism, it is difficult to real- 
ize an asymmetric particle accelerator in the turbulence 
region due t o frequent scatterings. Whether t he scenario 
proposed bv lMcClements fc Alexanderl (|2005l ) is the rule 
or an exception thus remains an open question. 

6. SUMMARY AND CONCLUSION 

We have presented imaging and spectral analysis of the 
RHESSI observations of the 2003 October 29 XIO flare 
showing two conjugate HXR footpoints (FPs), which are 
well-defined during the bulk of the fiare duration. One 
FP lies to the east (E-FP) and the other to the west (W- 
FP) of the north-south magnetic neutral line (NL). This 
fiare provides a unique opportunity to study in great de- 
tail the spatial, temporal, and spectral properties of the 
FPs and their associated magnetic fields. The impul- 
sive phase was relatively long (~20 minutes), HXR fiuxes 
were detected by RHESSI at energies up to hundreds 
of keV, and it was located close to disk center, result- 
ing in minimum projection effects and excellent magnetic 
field measurements from SOHO MDI. Our main findings 
regarding the unshearing motions, various correlations, 
and asymmetric characteristics of the two FPs are as fol- 
lows. 
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1. Two-phase FP unshearing and loop-top (LT) 
motions are observed in this flare. In Phase I the two 
identified FPs become closer to each other as they rapidly 
move almost anti-parallel to the magnetic NL, while in 
Phase II they move away from each other slowly, mainly 
perpendicular to the NL (Fig. [Th). In other words, the 
shear angle 6 between the normal to the NL and the line 
connecting the two FPs exhibits a fast and then slow 
decrease from 56° to 12° (Fig. [8^). This suggests that 
later reconnected magnetic field lines are less sheared 
(closer to a potential field), which is consistent with 
early observations in HXRs and other wavelengths (e.g., 
IZirin fc Tanakal [19731 ). More importantly, the transition 
between the two phases coincides with the direction re- 
versal of the apparent motion of the LT source along 
the NL (Fig. [3)), and the minima of the estimated loop 
length (Fig. Eh) and LT height (Fig. Eh). This suggests 
that the initial decrease of the L T altitude observed in 
many other RHESSI fiares (e.g., iSui Holm^^ [2003h 
may be associated with shorter loops during the fast un- 
shearing motion phase when the reconnection site prop- 
agates along the arcade. A possible explanation f or this 
early phase is the implosion conjecture (Hudsonl l2QQQl ) 
that predicts contraction of field lines during a solar ex- 
plosion including fiares and CMEs. 

2. There are correlations among the temporal evo- 
lutions of various quantities (Table [2]), some of which 
exhibit distinctions between the two HXR peaks (divi- 
sion at 20:48 UT): (a) The HXR fiuxes (Figs. [9]a and 
fTOb ) and spectral indexes (Figs. [9]c and [TUb ) of the two 
FPs are strongly correlated. This is evidence that the 
two HXR sources are from conjugate FPs at the two 
ends of the same magnetic loop, {h) The HXR fiux and 
spectral index of each FP show a commonly observed 
"soft-hard-soft" evolution (Figs. [9]a, [SH, and [TOfc ) dur- 
ing HXR Peak 1, while during Peak 2 the evolution be- 
comes "soft-hard-hard", (c) The magnetic field strengths 
at the two FPs also exhibit some temporal correlation 
(Figs.[9h andfTOb) particularly during Peak 1, consistent 
with the conjugate FPs identification, {d) The FP HXR 
fiuxes exponentially correlates with the magnetic field 
strengths fFigs. [TOle andfTOjf). which also anti-correlates 
with the spectral indexes during Peak 1. These correla- 
tions suggest that stronger magnetic fields, and/or larger 
reconnection rates or larger electric fields in the recon- 
nection region are responsible for producing larger fluxes 
and harder spectra for the accelerated electrons and thus 
the resulting HXRs. This is in qualitative agreement 
with the predictions of the stochastic acceleration model 
dPetrosian & Liu"2004') and the electric field acceleration 
model (Holman 1985). 

3. Various asymmetries are observed between the 
conjugate FPs (Table [3]): (a) On average, the eastern 
footpoint (E-FP) HXR fiux is 2.2 times higher than 
that of the western footpoint (W-FP; Fig. [Hh), while 
its magnetic field strength is 1.8 times weaker (520 G 
vs. 960 G; Fig. [9h). This is consistent with asymmet- 
ric magnetic mirroring ( ^5.ip . {h) The average estimated 
coronal column density from the edge of the LT source 
(assumed to be the acceleration region) to the transi- 
tion region at E-FP is 1.7 times smaller than that of 
W-FP (1.2 X 10^° vs. 2.1 X lO^^cm-^; Fig. [9]e). This 
qualitatively agrees with the HXR fiux asymmetry, be- 
cause a larger coronal column density results in more 



HXRs produced in the loop legs and thus less HXRs 
emitted from the FP below the transition region, espe- 
cially at low energies (Fig. [TTJ %.2\ . (c) The photon 
spectra above 50 keV of the two FPs are almost paral- 
lel to each other (Fig. [5]a), with the brighter E-FP be- 
ing consistently slightly harder than the dimmer W-FP 
(Fig. [9h). Their mean index values (71) = 3.63 ± 0.06 
and (72) = 3.79 ± 0.11 have a marginally significant dif- 
ference of (72) — (71) = 0.15 ± 0.13. In other words, 
the W-to-E ratio of the photon fiuxes is a constant or a 
slightly decreasing function of energy. This contradicts 
the column density effect which would produce a harder 
spectrum at the dimmer W-FP (Fig. [TTj). {d) As ex- 
pected from asymmetric magnetic mirroring, there is a 
temporal correlation between the W-to-E HXR fiux ra- 
tio and the E-to-W magnetic field ratio. However, this 
correlation only holds during HXR Peak 1 but gradually 
breaks down during Peak 2 (Figs.[9]c). We suggest that a 
combination of the asymmetric magnetic mirroring and 
column density effects could explain this variation ( ^5.3p . 
Specifically, since the column densities in later formed 
loops are larger (Fig. [9]e), collisions are more important 
at later times, making the HXR fiux ratio deviate from 
the prediction of mirroring alone. 

In our analysis we have treated the magnetic mirroring 
and column density effects separately in order to make 
the problem analytically tractable, yet without loss of 
the essential physics. However, in reality, the two effects 
are coupled and they should be studied together self- 
consistently to obtain a quantitative model prediction. 
This is do ne with the Fok ke r-Pla nck particle transport 
model of iLeach fc PetrosianI ()l98ll ) in a converging mag- 
netic field geometry. Results from such an analysis will 
be presented in a future publication. In addition to nu- 
merical modeling, we have started a statistical study of 
RHESSI fiares showing double FP sources that are close 
to disk center and thus have less projection effects. We 
hope to conduct future joint observations with RHESSI, 
Hinode, and the Solar Dynamic Observatory to obtain 
more advanced measurements of the magnetic fields at 
FPs. These future investigations will help improve our 
understanding of the underlying physics of asymmetric 
HXR FPs. 
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APPENDIX 

A. EFFECTS OF PULSE PILEUP ON IMAGING SPECTROSCOPY 

Pulse pileup refers t o the phenome non that two or more photons close in time are detected as one photon with their 
energies summed (Sm ith et al.ll2QQ2l ). When count rates are high, as happens in large flares, an artifact appears in 
the measured spectrum at twice or a higher multiple of the energy of the peak of the count rate spectrum that is at 
^6 keV in the i?iJE'S'6'/ attenuator AO state, ~10 keV in the Al state, and ~18 keV in the A3 state. Due to nonlinear 
complexity, there is currently no 100% reliable pileup correction algorithm available in the RHESSI software , especially 
for imaging spectroscopy. Here we have adopted and improved upon the methods used bv lLiu et al.l (|2006bl . their §2.1) 
to estimate the pileup importance and minimize its effects on our analysis. 

A general indicator of pileup severity is the livetime, the complement of the deadtime during which the detector is not 
able to distinguish among different incident photons. We find that the fractional livetime averaged over detectors 3-9 
and over 4 s intervals has a V-shaped time profile (not shown) with values >90% at the beginning (<20:40 UT) and end 
(>21:00 UT) of the flare, and a minimum of 24% at 20:46:10 UT during the impulsive peak. This value is very small 
compared with the livetime minima of 55% during the 2002 July 23 X4.8 flare and 94% during the 2002 February 20 
C7.5 flare, and indicates severe pileup effects. On much shorter time scales, the livetime fluctuates in anti-correlation 
with the count rate due to the modulation of the RHESSI grids during the spacecraft rotation, and a variation between 
5% and 70% is found for detector 9 at 20:46:10 UT. Such fine temporal variations make pileup correction for imaging 
spectroscopy even more difficult than for spatially integrated spectroscopy. 

ISchwarta (|2008) has recently developed a forward-modeling tool to simulate pileup effects for imaging spectroscopy, 
based on input model images consisting of sources of user-specified spectra and spatial distributions. This tool 
has been applied here to generate simulated images at 20:46:00-20:46:04 UT (near the livetime minimum) using 
the real flight instrument response. The model image consists of a thermal elliptic-gaussian LT source and two 
nonthermal (power-law) circular-gaussian FP sources, whose geometric and spectral parameters were selected to match 
the observation at 20:46:00-20:46:20 UT. We have run simulations for four cases as listed in Table [4] with different 
count rates corresponding to different levels of pileup. Case A returns the nominal result assuming no pileup at all; 
Case C represents this flare with the count rate of CR= 1.1 x 10^ (counts s~^ detector"^) identical to the maximum 
of the measured value at 20:46:10 UT. In each case images in 16 energy bins (used in §2.2p from 6 to 150 keV were 
reconstructed and then normalized as if their corresponding count rates were identical to the measured value. The 
final normalization makes different cases directly comparable and emphasizes the effects of pileup. 

TABLE 4 

Simulated effects of pulse pileup at different levels (indicated by the count rates) on imaging spectroscopy. The highest 

COUNT RATE DURING THIS FLARE IS CRobs = 1-1 X 10^ COUNTS DETECTOR"^. THE FLUXES AND EMS ARE normalized BY A FACTOR OF 

CRobs/CR, WHERE CR IS THE COUNT RATE ASSIGNED TO EACH CASE. 

Case Count Rate (CR) FP spectral index FP flux at 50 keV LT spectral parameters 

(cts s-i det-i) in 50-150 keV (ph s"^ cm"^ keV~^) in 6-20 keV 



(Subscript 1: E-FP, 2: W-FP) 


71 


72 


72 - 71 


h 


h 


h/h 


EM (10^9 cm-3) 


T (MK) 


A 7.7 X 10"^ (but no pileup) 


3.71 


3.82 


0.11 


25.5 


12.0 


0.47 


15.4 


30.8 


B 1.6 X 10^ 


3.74 


3.85 


0.11 


29.2 


13.7 


0.47 


15.0 


30.9 


C 1.1 X 10^ (measured) 


3.85 


3.93 


0.08 


52.8 


24.0 


0.46 


12.1 


31.0 


D 3.9 X 10^ 


3.66 


3.71 


0.05 


134 


58.9 


0.44 


6.64 


31.1 



Let us first examine the effects of pileup on source morphology. From Figure [12] we find the following features 
in agreement with the speculation of iLiu et"an (|2006bl . their §2.1): (1) The source positions, shapes, and sizes are 
generally well preserved, especially at low and high energies. This is because photons from the same region within the 
FWHM resolution of a grid have higher probability to pileup as one count in the underlying detector than photons 
from different regions. Therefore during image reconstruction piled- up photons are likely to be registered back to their 
original source. (2) The main visible artifact is the "ghost" LT source at intermediate energies (e.g., 34-41 keV). This 
is because the count-rate spectrum of this flare (in attenuator state A3) is ~18 keV at which the LT source dominates 
over the FP sources, and twice this energy is 36 keV where the first-order pileup is most pronounce. 

Next we check the effects of pileup on source spectra. In general, pileup shifts photons from low energies to high 
energies, as can be seen from Figures [T3l a and \Wh . For the LT source, the spectrum is only slightly reduced at low 
energies (^20 keV), with negligible changes in shape, while it is substantially altered in both normalization and shape 
at high energies. The opposite is true for the E-FP source whose spectrum >50 keV is only elevated by a roughly 
constant factor independent of energy (with minimal changes in slope). For comparison, we also perfo rmed pileup 
correction to the spatially integrated spectrum at 20:46:00-20:46:20 UT using the preliminary algorithm of'Smith e t al.l 
(2002) with a tweak factor of 0.6. The result shown in Figures [T3lc and [T3H. Because the LT source dominates at low 
energies, while the FPs dominate at high energies, the above mentioned features at these extreme energies are evident. 
The most obvious alteration of the spectrum shape occurs in the 30-40 keV range that covers twice the peak energy 
(18 keV) of the count-rate spectrum. 
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Fig. 12. — Simulated images for 20:46:00-20:46:04 UT at different energies, (a) and (6) Case A without pileup. (c) and {d) Case C with 
pileup representing this flare. The dark contours at 40% and 80% of the maximum in (6) and (d) represent the 6-8 keV image where the 
LT source dominates. 
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Fig. 13. Left: Simulated photon spectra for the LT (a) and E-FP sources (b) for Cases A, C, and D at 20:46:00-20:46:04 UT. The 
excessive LT flux >30 keV and E-FP flux <30 keV are due to the spatial overlap of the two sources. Right: Measured, spatial ly integrated , 
count rate spectrum (c) averaged over detectors 3-9 for interval 20:46:00-20:46:20 UT before and after pileup correction (Smith et^TESol). 
The broken lines are the estimated flrst- (two photons) and second-order (three photons) piled-up spectra, with their peaks marked by the 
vertical dotted lines, {d) The fraction of the total recorded counts that are piled up (including flrst- and second-order). The double arrows 
indicate the energy ranges used in the analysis. 

We repeated our imaging spectroscopy analysis here on the simulated data and the results are summarized in Table [H 
The striking feature is that the spectral shapes are well preserved across a wide range of pileup severity, as indicated 
by (1) the temperature T inferred from the LT spectra <20 keV and (2) the spectral indexes 71 and 72 and their 
difference, and the flux ratio h/h inferred from the FP spectra >50 keV. The FP portion of these spectral parameters 
are of our main interest in this paper. As to the absolute normalizations^ according to Case C, the estimated FP 
fluxes can be elevated by a factor of 2 and the EM can be reduced by 20% for this flare. These are the upper limits of 
changes as we used here the maximum observed count rate at 20:46:10 UT, and the changes are reduced considerably 
away from this time. The effects on the FP spectra are further reduced later in the flare as the spatial contamination 
from the LT is alleviated when the FP sources move away. We conclude that our choice of fltting the LT spectra only 
below 20 keV and the FP spectra only above 50 keV is justifled and has effectively minimized pileup effects on our 
analysis. We note in passing that the count spectrum of radiation belt particles estimated in § 12.11 appears as a hump 
peaking at ~50 keV and the counts integrated in the 3-600 keV band amount to only <8% of the total counts during 
20:40-21:05 UT. Therefore these particles have a negligible contribution to pileup. 

B. COALIGNMENT OF IMAGES OBTAINED BY DIFFERENT INSTRUMENTS 

We describe here how we c oaligned the RHESSI images and SOHO MDI magnetograms presented in this paper. 
(See a general tutorial: http://hesperia.gsfc.nasa.gov/~ptg/trace-align.) RHESSFs images are located on the Sun 
to sub-arcsecond a ccuracy thanks to its solar limb sensing aspect system and star based roll angle measurements 
(|Fivian et al.ll2002f ). MDI images had accurate plate scales and their roll angles were corrected for the solar P-angle, 
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but had different absolute origins for the solar x and y coordinates. MDI images were corrected to match the RHESSI 
features using the following procedures. 

(1) The first step is to identify specific features on the MDI map that have RHESSI counterparts. iQiu fc Garvl 
(|2QQ3l ) found good spatial agreement between HXR FPs and MDI magnetic anomaly features with an apparent sign 
reversal in a white-light flare. This was interpreted as HXR-producin^ nonthermal electrons being responsible for 
heating the lower atmosphere, which consequently altered the Ni I 6768 A line profile that is used by MDI to measure 
the magnetic field. 

We selected two neighboring magnetograms at 20:41:35 and 20:42:35 UT when the magnetic anomaly features 
were most pronounced, and subtracted the former from the latter. This running-differenced map (which we call 
mapo, Fig. [T4|) highlights regions of the newest changes, which are presumably caused by precipitation of nonthermal 
electrons and are expected to appear cospatial with HXR sources. As evident, there is one (three) apparent increase 
(decrease) feature(s) mainly in the negative (positive) polarity (cf.. Fig. [7j), which appear as white {dark) patches in 
Figure [HI Meanwhile, we reconstructed a RHESSI imdige (called mapi) at 60-120 keV integrated in the interval of 
20:42:19-20:42:51 UT (a multiple of the RHESSI spin period, ~4 s, and closest to the corresponding integration time 
of the second MDI magnetogram, 20:42:20-20:42:50 UT). We found an one-to-one correspondence between the four 
major HXR FP sources (Fig. [TH contours) and the magnetic anomalies. 

(2) We then converted SOHO^s LI view to the appearance as viewed from an Earth orbit, and we call the resulting 
differenced MDI magnetogram mapo, Earth- 

(3) Finally, we took the absolute value of MDI mapo, Earth to make a new map called mapQ ^arth coregister it with 
RHESSI mapi by cross-correlation. The required shifts for mapQ Earth — ^ — —2.8'' ± 2.0'', 
where the 2.0" uncertainty was estimated with the MDI pixel size. 



MDI diff. map, 20:42:35 - 20:41:35 U 
(x, y) shift= (4.5", -2.8") 



Fig. 14. — Difference {gray background) between two MDI magnetograms at 20:42:35 and 20:41:35 UT, overlaid with a RHESSI 60- 
120 keV image at 20:42:19-20:42:51 UT (white contours at 15%, 30%, & 80% of the maximum brightness). The MDI map has been 
corrected to an Earth-centered view and shifted in x and y to best match the four RHESSI HXR sources (see text). The dark dashed hue 
is the same simphfied magnetic NL shown in Fig. Oz. 



C. DERIVATION OF FOOTPOINT HXR FLUXES RESULTING FROM ASYMMETRIC CORONAL COLUMN DENSITIES 

Here we derive the numerical expressions for the HXR fluxes of the two FPs and their ratio as a function of energy 
resulting from column density asymmetry addressed in S I5.2[ We adopted the empirical expression of lLeach fc PetrosianI 
(| 19831 , their eq. 11) for nonthermal bremsstrahlung X-ray emission as a function of column density. This expression 
well approximates the Fokker-Planck calculation of particle transport under the influence of Coulomb collisions that 
includes energy losses and pitch-angle scattering, the latter of which was neglected in other forms of X-ray profiles based 
on approximate analytical solutions (e.g., Emslie & Machado 1987). For an injected power-law (index 5) electron flux, 
the resulting fractional bremsstrahlung emission intensity per unit dimensionless column density r at photon energy 
k (in units of rest electron energy mgC^ = 511 keV, mg being the electron mass) can be written as 

where r = A/" [47rrQ In A] = iV/[5 x 10^^ cm~^] is the dimensionless column density, for the classical electron radius 
ro = e^/rrieC^ = 2.8 x 10~^^cm and the Coulomb logarithm InA = 20; ^ is a factor determined by the pitch -angle 
distribution of the injected electron spectrum, which we assumed to be isotropic and thus g = 0.37 (|Leachlll983 ). This 
emission profile is normalized to unity, Io{r^k)dr = 1. Integrating equation ()Cip yields the cumulative photon 
emission from the injection site (r = 0) to the transition region (r = rtr = A^tr/[5 x 10^^ cm~^], where A^tr = Jq^"^ n[s]ds 
and Sir are the coronal column density and distance to the transition region), 

i^Corona(rtr) = ^0 (t, k)dr = 1 - (^1 + ^tr^J^ J , (C2) 
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whose complement gives the emission accumulated below the transition region, i.e., the HXR flux of the FP, 

roo / . -.x 1-5/2 

i^Fp(rtr) = J /o(r, k)dr = 1 - Fcorona(rtr) = (^1 + ^tr^ j • (C3) 

Note that at large photon energies (tens to hundreds of keV), Fcorona('rtr) is usually much smaller than FFp(Ttr). In 
addition, Fcorona(Ttr) is distributed in a large volume in the leg of the loop in the relatively tenuous plasma, while 
^Fp(Ttr) is concentrated at the FP in the dense transition region and chromosphere. This results in an even smaller 
surface brightness in the leg than at the FP, which may well exceed the dynamic idjU^e of HX R telescopes (e.g., ^10:1 
for RHESSI). This is why leg emission is so rarely observed (jLiu et al.ll2006bl : iSui et al.ll2006f ). 

As we know, a power-law electron flux (index=(5) produces a thick-target (i ntegrated from r = to r = oo) 
photon spectrum of approximately a power-law, /thick = Ao^~^, where 7 = ^ — 1 (Br ownll 19711 : lPetrosianl ll973) for an 
isotropically injected electron spectrum, and Aq is the normalization factor [in units of photons s~^ cm~^ (511 keV)~^]. 
Since /q gives the fractional spatial photon distribution at a given energy, the physical photon spectrum at energy k 
and at a depth where the overlying column density is r can be written as /(r, k) = /thick /o(t, k) = AoA:~^/o(r, k). It 
follows that the X-ray flux of the FP is 

roo / -r .x 1-5/2 

/Fp(rtr, k) = J I{r, k)dT = /thick/^Fp(rtr) = Aok-^ i^l + rtr^j , (C4) 

and the photon flux ratio of the two FPs (1 and 2), 

The above two equations were used in § 15.21 to calculate the FP fluxes and their ratio resulting from different coronal 
column densities. 

D. ESTIMATION OF COLUMN DENSITIES IN LOOP LEGS 

We describe below the approach to estimate the coronal column densities A/tr in the legs of the loop, which is defined 
and used in § 15.21 as the density integrated alon g the loop from t he ac celeration region to the transition region at 
the FPs. In the stochast ic acceleration niodel oflPetrosian fc Liul (|2004l ). the LT source is the region where particle 
acceleration takes place (|Liu et al.l 120081 : IXu et al" 20081 ). Assuming that this picture is true,^ we thus subtracted^ 
the estimated LT size (i.e., the radius r of the equivalent sphere; see Fig. [6]c) from the distances along the loop from 
the LT centroid to the FP centroids obtained in § 14.21 (i.e., where i=l for E-FP and 2 for W-FP; see Fig. [8]a), to 
obtain the path lengths in the legs Str,i = k — r. Here the FP centroids are assumed to be situated at negligibly small 
distances below the transition region. To give the desired column densities A/tr,i, the path lengths Str,i (^=1, 2) were 
then multiplied by the density nieg (assumed to be uniform) in the legs of the loop, which was estimated as follows. 

The density of the LT source tilt inferred in § 12.21 (see Fig. [6H) provides our first guess for the leg density nieg 
as assumed bv lFalewicz fc Siarkowski (2007). The relative brightness of nonthermal bremsstrahlung emission the leg 
and FP provides another important clue. This is because, for the same reason of collisional losses mentioned in § 15. 2[ 
the ratio of the leg to FP brightness, particularly at low energies, is an increasing function of the leg density. This 
predicted ratio cannot exceed the observed LT-to-FP brightness ratio, because the LT source is where the maximum 
loop brightness is located, and it includes additio nal contr ibutions from thermal emission, piled- up photons, and/or 
electrons trapped in the acceleration region (Petr osian fc Li u 2004). This imposes an upper limit for the leg density 
nieg, which we found to be nieg, max = 0.5nLT based on an error and trial method. This result indicates that the 
average density in the legs is smaller than the estimated LT density, or that the LT density is an overestimate due to 
an underestimate of the volume which could result from the lack of knowledge of the source size in the third dimension 
along the line of sight. This leg density was then used for column densities A/tr,i = ^ieg,max(^i — ^) shown in Figure [9]e. 
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